NASA TECHNICAL 


NASA TR R 432 



CLASSICAL SEVENTH-, SIXTH-, AND FIFTH-ORDER 
R IJ N GE-KUTT A-N Y ST R OM FORM IJLAS 
WITH STEPS! ZE CON FROL FOR GENERAL 
SECOND-ORDER DIFFERENTIAL EQUATIONS 


by Eric in fehlberg 

George C. Marshall Space flight Center 
Marshall Space flight Center , Ala. 35812 





NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. 


OCTOBER 1974 



1. REPORT NO. 

NASA TR R-432 


2. GOVERNMENT ACCESSION NO. 


3. RECIPIENT'S CATALOG NO, 


4. TITLE AND SUBTITLE 

Classical Seventh-, Sixth-, and Fifth-Order Range- Kutta-Ny strom 
Formulas With Stepsize Control for General Second-Order Differ- 
ential Equations 


5. REPORT DATE 

October 197*+ 


6. PERFORMING ORGANIZATION CODE 


7. AUTHOR{S) 

Erwin Fehlberg 


8. PERFORMING ORGANIZATION REPORT ft 

M546 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 


10. WORK UNIT NO. 


George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


1 1. CONTRACT OR GRANT NO. 


13, TYPE OF REPORT & PERIOD COVERED 


12. SPONSORING AGENCY NAME AND ADDRESS 


National Aeronautics and Space Administration 
Washington, D. C. 20546 


Technical Report 


14. SPONSORING AGENCY CODE 


15. SUPPLEMENTARY NOTES 

Prepared by Computation Laboratory, Science and Engineering 


16. ABSTRACT 


Runge- Kutta-Ny strom formulas of the seventh, sixth, and fifth order are derived for the 
general second-order ( vector) differential equation x = f( t, x, x) . These formulas include a 
stepsize control procedure, based on a complete coverage of the leading term of the local 
truncation error in x. The formulas require no more evaluations per step than our earlier 
Runge-Kutta formulas for x = f(t, x) . Since we have not to convert the second-order differential 
equations into twice as many first-order differential equations, the new formulas can be expected 
to be time-saving compared with the Runge-Kutta formulas for first-order differential equations. 
Two examples are presented. With results being of the same accuracy, in these examples the 
new Runge-Kutta-Nystrom formulas save from 25- to 60-percent computer time compared with 
our earlier Runge-Kutta formulas for first-order differential equations. 


17. KEY WORDS 


TTiT 


DISTRIBUTION STATEMENT 


CAT. 19 


19. SECURITY CLASSIF. (of thl» reports 

20. SECURITY CL ASSIF. (Of this page) 

21. NO. OF PAGES 

22. PRICE 

Unclassified 

Unclassified 

89 

$4.00 


MSFC - Form 3 2 9 2 (Rev December 1 9 7 2) For sale by National Technical Information Service, Springfield, Virginia 22151 




















ACKNOWLEDGMENT 


The author is indebted to Dr. Dale G. Bettis, Assistant Professor, 
Department of Aerospace Engineering and Engineering Mechanics, The 
University of Texas, Austin, Texas, for checking the equations of condition 
(Table l) by rederiving them in a different way and for applying the formulas of 
this report to various test problems. 

I also acknowledge the help of Mr. Barry Lee of the Computer Sciences 
Corporation ( CSC) , Huntsville, Alabama, for his programming work on our 
UNIVAC-1108 computer in the search for suitable coefficients for our formulas 
and for his programming in 40-digit arithmetic of the final coefficients (Tables 
3 and 5) as well as his computation in 30-digit arithmetic of the solution of 
Problems I and II. 

Finally, my thanks go to Dr. Helmut Hoelzer and Mr. Carl Prince, 
Directors of the Computation Laboratory, S&E, Marshall Space Flight Center, 
Huntsville, Alabama, for continuously supporting my work even under uneasy 
circumstances and for always making available to me the computing equipment 
of the laboratory. 



TABLE OF CONTENTS 


Page 


INTRODUCTION 1 

SECTION I. THE EQUATIONS OF CONDITIONS FOR RUNGE- 

KUTTA-NYSTROM FORMULAS 2 

SECTION II. SEVENTH-ORDER FORMULA RKN-G- 7(8) 9 

SECTION III. SIXTH-ORDER FORMULA RKN-G-6( 7) 27 

SECTION IV. FIFTH-ORDER FORMULA RKN-G- 5(6) 33 

SECTION V. APPLICATION TO TWO NUMERICAL EXAMPLES ... 39 

REFERENCES 41 


LI ST OF TABLES 


Table Title Page 

1. Equations of Condition for the Runge-Kutta-Nystrom 

Coefficients 42 

2. Pattern for RKN-G-7( 8)-13 ’70 

3. Coefficients for RKN-G-7( 8) -13 71 

4. Pattern for RKN-G-6( 7) -10 76 

5. Coefficients for RKN-G-6( 7) -10 77 

6. Pattern for RKN-G- 5(6) -8 80 

7. Coefficients for RKN-G- 5( 6) -8 81 

8. Application of the Various Formulas to Problem I 82 

9. Application of the Various Formulas to Problem II 83 


CLASSICAL SEVENTH-, SIXTH-, AND FIFTH-ORDER 
RUNGE-KUTTA-NYSTROM FORMULAS WITH 
STEPS IZE CONTROL FOR GENERAL SECOND- 
ORDER DIFFERENTIAL EQUATIONS 

INTRODUCTION 

1. In two earlier reports [1], [2] this author derived Runge-Kutta-Nystrom 
formulas for a special class of second-order (vector) differential 
equations 

x = f( t, x) ( 1) 

which do not contain the first derivative x on the right-hand side. In 
this report we will derive Runge-Kutta-Nystrom formulas for general 
second-order (vector) differential equations: 

x = f(t,x,x) (2) 

2. Similar to the Runge-Kutta-Nystrom formulas of report [ 1] and the 
Runge-Kutta formulas of report [3] , the formulas of this report include 
an automatic stepsize control based on a complete coverage of the leading 
term of the local truncation error in x. This coverage is achieved by one 
additional evaluation of the differential equations. Each of our Runge- 
Kutta-Nystrom formulas represents a pair of integration formulas for x 
which differ from one another by the one additional evaluation of the dif- 
ferential equations. The orders of these two formulas differ by 1. There- 
fore, the difference of the formulas represents an approximation of the 
leading term of the local truncation error in x for the lower-order 
formula. By requiring that this difference remain between preset limits, 
an automatic stepsize control for the lower-order formula can be 
established. 

3. The formulas for x are of the same order as the lower-order formulas 
for x. There is no automatic error control with respect to the formulas 
for x. Such a control would require a considerable additional effort, 
and such formulas would require more evaluations per integration step 
than our Runge-Kutta formulas of report [3] for first-order differential 
equations. 


However, when deriving the Runge-Kutta-Nystrom formulas of this report, 
we have also considered the error terms in x and we have selected the 
coefficients of our formulas in such a way as to keep the error terms in 
x as small as possible. 

SECTION I. THE EQUATIONS OF CONDITIONS FOR 
RUNGE-KUTTA-NYSTROM FORMULAS 


4. The derivation of the equations of condition for our Runge-Kutta-Nystrom 
formulas, as explained in this section, is based on a procedure of 

D. Sarafyan [4] . Sarafyan’ s method is extended to second-order dif- 
ferential equations to yield the equations of condition for the coefficients 
of Runge-Kutta-Nystrom formulas. 

5. In the following we explain in detail the procedure for a fourth-order 
Runge-Kutta-Nystrom formula. 

Let the evaluations for the Runge-Kutta-Nystrom formula for (2) be 


fn = 


f, = 


f, = 


f, = 


£(t 0 ,x 0 ,x 0 ) 

f ^t 0 + cqh, x 0 + x h + i f^ 2 h 2 , x 0 + 

f[to + « 2 h > x 0 + x 0 O! 2 h + + T 2 l( f l - f o) h2 > 

x 0 + f 0 Q! 2 h +/3 21 (fj - f 0 ) h J 

f jt 0 + a 3 h, x 0 + x 0 a 3 h +| f 0 o^ h 2 + y 31 (f! - f 0 )h 2 + y 32 (f 2 - f 0 )h 2 , 
x 0 + foc^h + /3 3 i(fi - f 0 )h + /3 22 (f 2 - £o)hJ 


(3) 


2 


I 



with t 0 ,x 0 ,x 0 being the initial conditions for the integration step, h the 
stepsize and the a's, /3’s and y's the Runge-Kutta-Nystrom coefficients 
that we want to find. 

The evaluations (3) lead to the Runge-Kutta-Nystrom formulas: 


x(t 0 + h) = x 0 + x 0 h + (c 0 f 0 + c 4 fj + c 2 f2 + c 3 f 3 + c 4 f 4 + c 5 f 5 )h 2 + 0(h e ) 
x(t 0 + h) = x 0 + (c 0 f 0 + Cjfj + c 2 f, + c 3 f 3 + c 4 f 4 + c 5 f 5 )h + 0(h 5 ) 


(4) 


with weight factors c and c that we also want to find. 

K K 

6. To determine these unknown coefficients, we expand the solution x(t) of 
(2) into a Taylor series at t = t 0 : 


x(t) = x 0 + x 0 (t - t 0 ) + A 2 (t - t 0 ) 2 + A 3 (t - t 0 ) 3 + A 4 (t - t 0 ) 4 
+ A 5 (t-t 0 ) 5 + O(t-t 0 ) 6 


( 5 ) 


with the abbreviation: 


A = 

K 


k! 


djc\ 

dt K / 0 


_1 

k'. 



(« 


2,3, . . .) 


( 6 ) 


Inserting (5) into (2) and setting t - t 0 = h yields: 

f(j + 6A 3 h + 12A 4 h 2 + 20A 5 h 3 + . . , 

= f(t 0 + h, x 0 + x 0 h + A 2 h 2 + A 3 h 3 + A 4 h 4 + A 5 h 5 + . . . , 

x 0 + 2A 2 h + 3A 3 h 2 + 4A 4 h 3 + 5A 5 h 4 + . . . ) = R 


( 7 ) 


For the sake of briefness, we denoted the right-hand side of (7) by R. 
Taylor expansion of R at 


(0) : (t 0 + h, x 0 + x 0 h, x 0 + 2A 2 h) 


( 8 ) 
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leads to 


R = 




( A 2 tf+A 3 tf ) 

^3A o h 2 + 4A 4 h 3 ^ + 0(h 4 ) 


(9) 


An expansion at 

0 : (t 0 ,x 0 ,x 0 ) (10) 


leads to 



introducing the operator: 

' D <». =[FT() + ^O ^^ ()t ] 0 • (12) 


We insert (9) and (ll) into (7) and replace h by a h: 

K 



4 



Because of (4) the functions f in (3) must be correct up to the third 
power in h. We set 



with constants q„ , q„ , p„ , p, that will be determined as follows 
m 0k’ Ik Ok Ik 

somewhat later. 

By Taylor expansion of f at 

K 

<0> : ( to + a h, x 0 + x 0 a K h, x 0 + f^li) (15) 


we obtain 


f = 


(f ) 

' k' 


< 0 > 


/of \ . 

+ ( — — ) / q h 2 + q h 3 ) 

\9 x/ 1 4 0 k 4 1k ) 

' ' < 0 > \ / 

© / V 

( p h 2 + p h 3 ) 

^ V 0K lK / 


(16) 


We now insert (13) and (ll) into (16) after having replaced h in (ll) 

by a h: 

J K 


f = f + 6A a h + 12A a 2 h 2 + 20ha 3 h 3 
k 0 3k 4k 5k 


+ (q - A a 2 \h 2 + (■£—) (q, - A a Ah 3 

\9xy o y Ok 2 Kf yaxy^y Ik 3 kJ 

+ [ D (^)] 0 “'<( q » K " A 2“«) hJ 


(17) 


5 


(U ) 0 *(H) ( p i„ ■ 


[Ksi 


<>(P^ - 3A 3 a 2 )h 3 


Ok 


+ 0(h 4 ) 


(17) 
(con. ) 


Equation (17) is called the generating formula since the introduction of 
(17) into (4) generates the equations of condition for the Runge-Kutta- 
Nystrom coefficients. 

7. For k = 1 we find from (14) and the second equation (3) 


1 2 

Qoi = g f o«i » Qii = 0 . Pol = 0 . P 11 = 0 


(18) 


and from ( 17) 

f t = f 0 + 6A 3 0! 1 h + 12A 4 a! 2 1 h 2 + 20A 5 ct!jh 3 - (■^—) A 3 o; 3 h 3 

\ax/ 0 

- 3 (H ) 0 A » tt > h! - 4 (H)„ A <“- h3 - 3 [ D (H)], A »“' h3 

+ 0(h 4 ) 


(19) 


For k = 2 we find 


Q _ 7T fo a 2» Ql2 " 6A 3 y2iQ!i 
02 z 


/ 0 A 

P02 = ®A 3 ^2l 0! l> Pl2 — 12A^21^1 - g x/q A 3 ^21 Q! l 


( 20 ) 


We introduce (20) into (17) and omit all terms that have already corre- 
sponding terms in ( 19) : 
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u = 


• + 6 (H) o A 3T2l Q! l h3 + 6 (H) 0 A 3^2l Qi l h2 

+ 12 (H ) 0 A A\ a \ h3 - 3 (||) o A 3^h 3 

D (fi)] 0 A ^2i a l h3 + °( h4 ) 


+ 6 


Continuing, we find for k = 3: 


( 21 ) 


Q03 = | f 0 a 3 » QlS = G A 3(Yn a i + T 3 2« 2 ) » P03 = + 032“*) 

/ 0 £\ / 3 f. \ 

Pl3 = 12A 4 (/3 31 0f;j + /^32 ce 2 ) - 3 ( g"x ) 0 ^(^l 01 ! + P'A&l) + ) A !1^32^21 G! 1 


( 22 ) 


and again omitting all terms that have already corresponding terms in 
(19) or (21): 


( b f\ 2 
— . . . + 61 — r I 

\3x/ ( 


A 3 / 3 32/ 3 2i a l h3 + 0(h 4 ) 


(23) 


No more additional terms are obtained for k = 4 and k = 5. 


8 . 


The equations of condition for the Runge-Kutta-Nystrom coefficients are 
obtained by inserting the expansions (5) and ( 19), (21), (23), ... into 
the Runge-Kutta-Nystrom formulas (4). 


From the first equation ( 4 ) we find 
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*o + x 0 h + A 2 h 2 + A 3 h 3 + A 4 h 4 + A 5 h 5 

x 0 + x„h + (c 0 + . . . + e s )f„h 2 + 6A 3 (C|Qfj + . . . + CgG^h 3 

12A 4 (c 1 oi 2 i + . . . + c 5 a|) h 4 + 20A 6 (c,<* 3 + . . . + c 5 af)h 5 


at 


j£). + • • • + c 5«l)h 5 - 3(||) A 3 (c,a? + . .. + c s Of|)h‘ 


' 4 \8x/ 0 A4(c >“ ? + ' • + c 5“i )h5 - 3 Pin?) L + • • • + c s <*?)h 3 

/ g j\ 

+ 6 (ax j 0 A 3 [c 2 y 21 Q! i + .•. + Cgfy^Q!, +y 52 o ; 2 + 753^3 + y 54 0 ! 4 ) } • h 5 

+ 6 (H) 0 A lI«Al“l + --- + c,(0,,a, + 0^ 2 +0 M a 4 )] • h 4 

/ 8 f \ 

+ 12 (axy 0 A4fc * -M*..- + 0 5 (/3 5I o5 + /J 52a 2 + /3 53 « 2 +/ 3 54a ,2)] h 5 

/ ^ f \ 2 

“ 3 \9x/o A;j[C2 ’ + * * * + c *(P*l a \ + ^r , 2 «2 + Pm<xI + ^ 54 ^ 4 ) 1 h 5 


+ 6 P(57;L A3[C * af2 ' + C S^ 5 1*1+/V*2 +^3 + ^54«4)ih 5 


( 9 f 

8xJ„ A 3 {c 3 ■ 0 32 . 0 2i or, + . . . + C 5 [3 52 . /3 21 a, + P 53 (/3 31 Q., + £J 32 o; 2 ) + /3 54 (0 41 a, + /3 42Q , 2 + /} a 


)l}h 5 


~\ 


( 24 ) 


J 


Equating corresponding terms in (24) yields 
C 0 + C 1 + c 2 + C 3 + c 4 + c 5 = i 

C 1 A 1 + 02^2 + C3O3 + C 4 0 ; 4 + C 5 Ct 5 = — 

b 

C l a \ 4 c 2 a 2 + c 3 a 3 4 c 4 a?l + c 5 o;i = ~ 
c i a l + + c 3 o| + c 4 a 3 4 + c 5 a| = ± 


(25) 


and 


c 2 Y 2 io;i + c 3 (y 31 of I + y 3a a 2 ) " e 4 (y 4l ot l * y i2 a z + y 43 « 3 ) + c 5 (v 5 ia t + y S2 a, + y 53 a' 3 + y 54 a 4 ) = —■ 

c 2^21 Q! i + Cs(^3J a J 4 ^32 rt 2) 4 c 4(^ll a l 4 ^42 a 2 + 4 c &(^!j! q; 1 4 ^52 a 2 ' r ^S3 a 3 ^ >*54 a l) = 7^ 

*hPu<*\ 4 c 3 (0 31 a? + *? 32 «| ) + c 4 (/? 4l or? +0 n a 2 +Pa a D 4 c ^hi a \ 4 P n a l 4 PaPl 4 ^54 a 5) = ^ 

c 2 o 2 • P 2 i a 'i 4 C3a 3 (^3i«i + i* 32 a 2 ) ■» c 4 or 4 (/ 3 4 |re l + 0 42 a 2 4 ^43^3) 4 c 5 a 5 (P s iari * /3 5 j°2 4 4 = ^ 

C 3^32^21 a ’l 4 C 4^42 ’ >’21 a l + »*43(^3t ar l 4 ^32 Ct 2 ) ^ 

4 c !jl^52 ' ^21 a 't ■* f^53(^31 a ' I 4 ^32°^) T »^j(^4l a I 4 P 12 Q 2 4 ^43^3) ^ ~ 7777 


( 26 ) 
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Similar equations of condition are obtained from the second equation (4). 
Equations (25) and (26) are listed as the first nine equations of Table l. 1 
The left-hand part of Table 1 represents the equations of condition for x; 
the right-hand part represents those for x . For the right-hand part the 
weight factors c have to be replaced by c as indicated at the top of 
the table. 


9. The procedure described in this section can be extended to cover higher- 
order terms in the Taylor expansions. The resulting equations of con- 
dition up to the ninth order for x and the eighth order for x are listed 
in Table 1. Their derivation is naturally somewhat more involved than 
in the case of a fourth-order formula, but it follows along the same lines 
and is rather straightforward. 

To shorten the equations we introduced in Table 1, the abbreviations 


« A A 

|3 a + /3 a + 

k1 1 k2 2 


X A 

y a + y a + 
y K 1 1 y k2 2 


+ 13 a = P . 

K,K-1 K-l kA 


+ y a A = Q 
y k , k - 1 K-l kA 


(27) 


SECTION II. SEVENTH-ORDER FORMULA RKN-G-7(8) 


2 


10. We shall present in the following a seventh-order formula based on 

thirteen evaluations, a fourteenth evaluation being taken over as first 
evaluation for the next step. 

Let the evaluations be 
fo = f (to» x o»^o) 

( ' K=1 K-l 

t 0 + a h, x 0 + x 0 a h + h 2 . £ y • f^, x 0 + h . £ 

A=0 A=0 

(k = 1,2,3, ... , 13) 


1. All tables are at the end of this report. 

2. We insert in the names of the formulas of this report the letter G, to indicate 
that these formulas hold for the general differential equation (2) in contrast to 
the formulas without G of [1] and [2] , which hold for the special differential 
equation ( l) . 



9 



and the Runge-Kutta-Nystrom formulas 


12 


x = x 0 + x 0 h + h 2 • Yi c K f K + °( h8 ) 


K = 0 


13 

X = X 0 + x 0 h + h 2 • Yj ^ K f K + °( h9 ) 


12 

x = x 0 + h- Y c f + 0(h 8 ) 

rV rv 


( 29 ) 


The first formula (29) is a seventh -order formula for x ; the second 
formula is an eighth-order formula for x. The difference x - & will 
represent a first approximation of the local truncation error for x and 
will be used as stepsize control. The third formula (29) is a seventh- 
order formula for x . 

11. Similar to our previous reports [1], [2], [3], we make a number of 
assumptions for the Runge-Kutta-Ny strom coefficients that will reduce 
the number of equations of Table 1 to such an extent that we can handle 
the remaining problem with relative ease. 

Let us assume 


Cl = 

Cj 0 j • * • 

> c 6 — c 6 

= 

0, <3 7 = c 7 , 

• • • 

» ^11 

c il> c i2 

0, 

Cl = 

0, • • • , c 6 

= 0 ; a 

12 

= a 13= 1 





*31 

= ^41 = * ’ ’ 

= /3 131 = 

0 

Y 31 = 

y 41 = 

= . . . 

- y =0 

r 131 


^52 

= ^62 = " * " 

= /3 132 = 

0 

Y 62 ~ 

Y 72 = 

= . . . 

- y - 0 

r 132 


^73 

“ ^83 “ ’ “ 

_/3 133 = 

0 

^13, A 

= c 

A 

(a = 

0,1,2, 

, 12) 

^74 

= ^84 = • - 

= ^134 = 

0 

Y 13,A 

= °A 

(A = 

0,1,1, ... 

, 12) 


(30) 


(31) 
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The assumptions (30) mean that only the last two weight factors of the 
first two formulas (29 ) differ. The last two weight factors in these two 
formulas are simply exchanged: c 12 , 0 is replaced by 0, c 12 . There- 
fore, we have to compute only one set of weight factors for these two 
formulas (29). 

The assumptions (31) are necessary in connection with the assumptions 
of No. 12 and No. 13 to reduce the equations of Table 1. 

12. Let us further assume 


1 

^ ~~ 2* K “ 2, 3, ... , 13) . (32) 

X— 1 


As one can easily verify, assumptions (32) eliminate a large number of 
equations of Table 1 by converting them into other equations of this table. 

In the following we list the equations of Table 1 that are eliminated by the 
assumptions (32) 

IV :2 
V:2, 5 

VI:2, 7, 8, 9, 10, 13 

VII:2, 7, 9, 11, 12, 13, 14, 17, 20, 25, 26, 27, 28, 29, 30, 31, 34 
Vin:2, 7, 9, 11,14,17,19,20,21,22,25,28,33,34,35,40,41,46,48, 
50,51,52,53,54,55,56,57,58,59,60,61,62,65,66,67,70,71, 
72,75,78,83,84,85,86,88,90,91,92 

Since the ninth-order equations (IX) of Table 1 enter the computation only 
as eighth-order equations for x , we will consider these evaluations when 
dealing with the local truncation error terms in x . 

We next assume 


K -1 

l 


X=1 


> a* 
kX X 



(k = 2,3, ... , 13) 


(33) 
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thereby eliminating the following equations of Table 1: 


V :4 

VI :4, 12 

VII:4, 16, 19,22,33 
VIII:4, 15, 18, 24, 27, 30, 37,43 , 
64,69,74,77,80,94 

The assumptions 


K-l 

z 


A=1 


ft o;^ 
kA A 



( K = 5,6, ... ,13) 


(34) 


eliminate the following equations: 
VI :6 

VII:6, 24 
VIE :6, 32,45 

and the assumptions 


i (K=7 - 8 13) (35) 


the equations: 

VII:10 
VIII: 10 

Finally we assume 

c 7 /3 75 + c 8 /3 85 + c 9 /? 95 + c 10 /3 105 + + c 12 /3i 35 = 0 

C 7 /3 T5 + C 8 0 85 + C 9 ^ 95 + C 10 /3 105 + Cu/lus + C 12 /3 125 = 0 


(36) 
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and 


c 7^76 + C 8^8G + c 9^96 + c 10^106 + c 11^116 + c 12^136 — 0 j 

C 7 ^ 7 S + c 8 /3 8 g + c 9 /3 96 + c 10 /3 10 g + c n /3 lle + c 12 /3 l2 g =0 j , 

thereby eliminating equations VIII:49, 82, 95 from Table 1. 
13. We make similar assumptions for the coefficients y : 


K_1 

Z 


A=1 


y ~ oj 

kA A 



(k = 2,3, ... , 13) 


(37) 


(38) 


eliminating 
V :3 

VI: 3, 11 

VII:3, 15, 18,21, 32 
Vin:3, 13,23,26,29,36,42, 
63,68,73,76,79,93 ; 

K-i 

Y. y, a l=77r a4 («= 5,6, ... , 13) 
^ - kA A 12 K v ' 

A— 1 


(39) 


eliminating 
VI: 5 

VII:5, 23 
VIII:5, 31,44 ; 


k- 1 

Z 


A=1 


y » cv? 
K A A 



( K = 5,6, ... , 13) 


(40) 


eliminating 
VII: 8 

VIII:8, 47, 89 


13 



and 


K -1 

E 


\=1 


y 

k X X 


30 


cr 


K 


(*=7,8, ... , 13) 


(41) 


that eliminates VIII:12 from Table 1. 
Finally we assume 


c 7'T / 73 + c 8T83 + c 9l / 93 + c loT / 103 + c ll7ll3 + c 127l33 ® 

C 7 T 73 + C 8 783 + ^9793 + Ci 0 7103 + Cil7il3 + C 12 7l23 = 0 


(42) 


C 7 7 74 + C 8 7 84 + c 97 94 + c lo7 104 + c ll>Tl4 + c 127l34 _ 0 
c 7 y 74 + c 8 y 84 + c 9 y 94 + c 10 y 104 + c 41 y 114 + c 12 y 124 ~ 0 


(43) 


thereby eliminating VIII:38, 39, 81, 87 from Table 1. 

14. The assumptions of No. 11, 12, and 13 reduce the eighth- and lower-order 
equations of Table 1 to the following equations: 

(11,1), (III, 1) , (IV, 1), (v, i) , (VI, 1), (VII, 1), (vm,i), (VIII, 16) 

These equations have now to be solved together with the assumptions (30) 
through (43). 


From equations 

(III, 1) c 7 a 7 + c 8 a 8 + c 9 g! 9 + c 10 a! 10 + + c 12 = - 

(IV, 1) C 7 Oi| + CjjG'jj + C 9 0! 2 9 + C 10 G! 2 10 + CjjQ! 2 !! + c 12 = 

(V, 1) c 7 o: 7 + c 8 « 8 + c 9 o; 9 + c l0 Q!f 0 + c^o; 3 !! + c i2 = 


(44) 


14 



(44) 
(con. ) 


(VI, 1) c 7 af + c 8 a 8 4 + ego?! + c 10 af 0 + c u a u + C 12 = ^ 

(VII, l) c 7 a Y + c 8 a 8 + CgO?! + c 10 a| 0 + c 11 o'| 1 + C 12 = 

(VIII, 1) c 7 af + c 8 o; 8 + CgQi | + c 10 a^ + cnQ! 6 ^ + c 12 = 


we find the weight factors c 7 , c 8 , c 9 , c 10 , c u , c 12 as functions of 
®8> ®9> a 10» ®11* 


From the right-hand sides of Table 1 we obtain for the corresponding 
equations 


(III, l) * C 7 a 7 + CgGJjj + Cgttg + C 10 a 10 + 6^0! tl + c 12 = - 

(IV, 1) * c 7 a 2 7 + c 8 o:| + Cg^i + c 10 a 2 10 + c^oi 2 !! + c l2 = | 

( V, 1) * c 7 a? + c 8 0 ! 8 3 + c 9 ag + c 10 a?o + 0 ^ 0 !^ + c l2 = ^ 
(VI, 1) * c 7 a 7 + Cgag 4 + CgO!g + c 10 oif 0 + c n a A n + c 12 = | 

(VII, 1) * c 7 a| + c 8 a 8 5 + CgO!g + c 10 af 0 + cn^fi + c 12 = ^ 


(VIII, 1)* c 7 af + 


C 8 «8 


CgOil + C 10 <^0 + CuQ! 6 !! + C 12 = - 


(45) 


the weight factors c 7 , c 8 , c 9 , c 10 , c^, c 12 as functions of a 7 , ce 8 , oig, 
a 10 , a n . Equation (II, l) yields c 0 , or c 0 when written as (II, l) * . 

15. We still have to satisfy equation (VIII, 16) and the assumptions of No. 12 
and No. 13. 

From (32) and (33) we obtain 
K — Z K — Z 



(46) 
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and as restrictive condition 


a i = 



( 47 ) 


In the same way we obtain from (32) 


k=3 


and (33) k=3 


^32 



(48) 


and 


a 


2 



Equations (32) 

K = 4 


and (33) yield 


^42 


1 , 3a o - 2a, 

— a^ — -V 

6 a 2 (a 3 - a 2 ) 


(49) 


(50) 


and a corresponding formula for /3 43 , obtained from 0 42 by exchanging 
a 2 and a 3 . 

Equations (32) , (33) , (34) yield 

rC — D rC — 0 K — O 


£53 


1 2 30!, - 2a c 

— ai — -A \ 

6 a 3 (a 4 - a 3 ) 


(51) 


a corresponding formula for j3 54 and the restrictive condition 


1 4a, - 3a s 

2 5 3a 4 - 2a 5 


Equations (32) , (33) , (34) give 

K-b K-b K=b 


^63 


_L a 2 6a 4 a 5 - 4(a 4 + a R )a fi + 3a| 
12 6 a 3 (a 4 - " 3 ) (“5 - “ 3 ) 


and corresponding formulas for j 3 S4 and /3 65 . 


(52) 


( 53 ) 
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Equations (32) (33)^_ 7> (34 ) and (35) represent four linear 

equations for the two coefficients /? 76 , /3 76 which lead to 


o 1 9 3 etc - 2 q; 7 

Bn, = — Oti t -2 4 - 

75 6 aJo^Za 5 ) 


(54) 


a corresponding formula for jS 76 and to the two restrictive conditions 


“ 5 = > a e = Jq ( 6 +^" 6 ) 


ce, 


(55) 


Equations (32) K=g> (33) (c=q , (34) K _ 8 * and ( 35 ) K=8 lead to 


„ 1 -.2 6o; fi Q; 7 - 4(a fi + a 7 )o! 8 + 3a? 

* 8B ~ 12 V(« 6 - a B V(a T - a 5 ) ^ ’ 


(56) 


corresponding formulas for /3 86 , /3 87 and to the restrictive condition 


a 8 = 4 ff ? 


(57) 


^95 


The four equations (32) k = 9> (33) K=9 » ( 34 ) K = 9 > and ( 35 ) K = 9 yield 

1 2 30Q! fi Qi 7 Q; R - 20(a; fi Q! 7 + ctfiCtj) + q? 7 a? a )o; 9 + 15(o; fi + a; 7 + a? R )Q; 9 - 120;! 

60 “ 9 « 5 (“6 -“5)(“7 - Q! 5)(“8 -“ 5 ) ~ * 

(58) 


and corresponding formulas for /3 9g , 0 97 , and (3 98 . 

Putting 

ft 105 = 0 , (59) 
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we obtain from (32) K = 1() , ( 33 ) K = 10 » ( 34 ) k = io’ and ^ 35 ^« = 10 : 

1 2 30q? 7 o' s q? 9 - 20(a! 7 Q' fl + g 7 0f 9 + + 15(a 7 + a fl + QiflMn - jjg|n 

Pl °e 60 ai ° a«( a 7 - “e)^ - a e)( a 9 - “e) 

(60) 

and corresponding formulas for 0 1O7 , /3 108 , /3 l09 . 

From the first equation (36) and the first equation (37) we obtain 


0115 ~ - ~ ( c 7075 + c 8085 + c 9095) 

C 11 

0116 = — ( c 70 76 + c 8086 + c 9096 + c 10 0 106 ) 

C 11 

and from equations ( 32 ) ( 33 ) ^ 34 ^k = ii’ and ^ 3 ^^« = li : 


(61) 


B = J_ „* SOttgttaam - 2^a t a t + inKi + 15(«» + «» - 12a ll 

p<1? 60 11 a 7 (a 8 -a 7 )(a 9 -a 7 )(a 10 - a 7 ) 


- 0115 


-Q!!i)(«a -«s)(«m - <*s) 
a 7 (ci ! 8 -a 7 )(a 9 -Q! 7 )(q ; 10 - a 7 ) 


3 . <*«(<*« -«fi)(Q!» -« B )(ain -«b) 

116 a 7 (a 8 -a 7 )(a 9 -a 7 )(a 10 - a 7 ) 


(62) 


and corresponding formulas for /3 118 , /3 119 , and 0 111O . 

Equations (32 ) k = 11 , ( 33 ) K=11 > ( 34 ) K= n’ ( 35 ) K=11 » and(VIII,16) can 

be considered as five linear equations for the four coefficients 0 117 , /3 j 18 , 
/3 119 , and /3 1110 . Therefore, a restrictive condition for the a’s can be 
derived from these five linear equations: 


a li 


N(ttii) 

D(«ll) 


( 63 ) 
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with 


N(an) = 70o' 7 a 8 Q! 9 0'io - 42(o; 7 o! 8 o! 9 + + o' 7 a 9 o' 1(l + a-ga^io) 

+ 28(o! 7 a 8 + d 7 d 9 + Q! 7 O! 10 + a 8 Q! 9 + ag^io + “^lo) 

- 20(a 7 + a 8 + a 9 + a 10 ) + 15 

- 840c 12 ( 1 -a 7 )(l -« 8 )( 1 - ® 9 ) ( 1 -«io) 

D(an) = 14OG' 7 G' 8 G' 9 a! 10 - 70(a' 7 o' 8 o' 9 + a 7 oi g o; 10 + o; 7 O! 9 o; 10 + Q! 8 O! 9 Q! 10 ) 
+ 42(a 7 G! 8 + 0! 7 q; 9 + g! 7 O! 10 + a 8 Q! 9 + o^io + “s^io) 

_ 28(ce 7 + 0! 8 + a 9 + a 10 ) + 20 

- 840c 12 ( 1 -a 7 )(l -« 8 )( 1 -“eK 1 - a io) 


and 


c 12 “ 


1 

28 


N( c tt>) 

0(012) 


(64) 


with 

N(c 12 ) = 70 Q! 7 Q! 8 0! 9 Q! 1o - 28(a 7 0' 8 Q! 9 + Q! 7 O; 8 Q' 10 + Oi id apt 10 + O^^io) 
+ 14 (a 7 a! 8 + d 7 a 9 + d 7 d lQ + d 9 d 9 + a 8 a 10 + a 9 a; l0 ) 

- 8(gi 7 + a 8 + o; 9 + a 10 ) + 5 

D(c 12 ) = 30q! 7 O! 8 o; 9 g! 10 - 20(a 7 a s a 9 + d 7 d s a i0 + d 7 d 9 a 10 + Q! 8 Q! 9 O! 10 ) 
+ 15(o; 7 a 8 + d 7 d 9 + a 7 a 10 + a 8 a 9 + ^“io + “^io) 

- 12(a 7 + oi 8 + a 9 + o; 10 ) + 10 
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The second equations (36) and (37) yield 


^125 


- (c 7 ^75 + C 8 ^ 85 + 69^95 + 

c 12 


0126 


- -r— (c 7 /3 76 + c 8 /3 86 + 69/399 + c 10 /3 106 + c 17 ^ 116 ) 
c 12 


( 65 ) 


We then obtain from (32) 
(VIII, 16)': K = 12 


(33) 


K = 12 ’ 


(34) 


K = 12* 


(35) k =i 2 , and 


0127 


n(£m) 

a 7 («8 - « 7 ) (<*9 - “7) (“10 - ® 7) (ail - a?) 


( 66 ) 


with 

N ( 0 12 7 ) = g^- [30 o' 8 Q'9Q! 10 q; 11 - 20(a 8 a 9 o( 10 + a & a 9 a u +ar 8 a 10 a 11 
+ o^io^n) + 15(o! 8 a9 + oi 8 a 10 + a^an + a 9 a! 10 
+ CKgOi^ + O'jqQ'h) - 12(ci ! 8 + 0 ( 9 + OL 10 + 01 lt ) + 10] 

- 0 125 Q! 5( Qf 8 - a s) ( a 9 - “5) ( a 10 * a s) ( a ll - a 5) 

- 0i26 a 6( a 8 - a e(( a 9 - a e)( Q! io -®e)( a n - “e) - R 

and corresponding formulas for /3 128 , /3 129 , /3 1210 , and /S 1211 . 

The abbreviation R stands for 

- (0125^5 + 0 126^ 6 ) - ~ ^ - C8087 q: 7 " ^9(097^7 + 098^1) “ ^10(0107*^7 

K L I 

0108*^8 + 0109*^1) “ Cn(0ii 7 « 7 + 0118*^8 + 0119*^9 + 0111O CV lo) 

( 67 ) 

This concludes the computation of the coefficients /3^. 



16. The computation of the coefficients y ' proceeds in a similar way. 

From (38) _ we obtain 

K Z 

y 2 i = \ » (68) 

and from (38) : 

K — o 

y 3 2 = \ a i • (69) 


I 


Putting 


T42 = 0 


(70) 


equation (38) , yields: 

k= 4 


T43 


I . 

6 <*3 


Equations (38) (39) k=5 , (40) k = 5 yield 


y 3 2 


J_ 10O! 3<*4 ~ 5(o?3 + <Xi)Qls + 30 
60 5 3 - a 2 )(Q! 4 - a 2 ) 


(71) 


(72) 


and corresponding formulas for y 53 , y 54 . 

In the same way, we obtain from (38) k _ 6> (39)^ , and(40) K _ 6 : 


^63 


J_ ^3 10^4^5 - 5(tt 4 + tt B )Q!6 + 3®l 

60 6 « 3 («4 -° ; 3)( a 5 " <* 3 ) 


(73) 


and corresponding formulas for y g4 , y 65 , and from (38) , (39) 

< 40 >«,7- a " d (4 I >k=7 : “ " 
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TT3 = 6? 


10ct 4 q? 5 q: 6 - 5(o' 4 Q' 5 + ce^ + Q' 5 0! 6 )o' 7 + 3 (ce 4 + o? 5 + Q' fi )o'f - 2a f 
«' 3 ( Q! 4 - " 3 ) («5 - «s) ( a 6 - « 3 ) 


( 74 ) 


and corresponding formulas for y 74 , y 75 , y 76 . 
Putting 


783 = 0 


(75) 


^84 


equations (38) K=g , (39) K==g , (40) K=g , (41) ^ =g yield 

J_ ^3 IQo^fiO'T - 5 (q' s Q'( ! + Q! S 0 ' 7 + Q' fi o; 7 )o; 8 + 3(a s + + Q' 7 )«^ - 20$ 

60 8 a 4 (« 5 - °U) (“e” - ^ 4 ) ( a 7 -<* 4 ) 


(76) 


and corresponding formulas for y 85 , y 86> andy 87 . 
Putting 


T 93 = T 94 = 0 

equations (38) k=9> (39) K=g , (40) K=g , and(4l) K=g lead to 


(77) 



10 o' 6 o! 7 a 8 - 5(o' fi a' 7 + Q! 8 o; 8 + q' 7 0 ! 8 )o' 9 + 3(q? 6 + aj + o' 8 )q; 9 - 20?! 
a 5 (Q!6 - « 5 ) («t " «s) ( a 8 - “ 5 ) 


(78) 


and corresponding formulas for y 9e , y 97 , andy 98 . 
With 


T 103 ~ y\04 ~ ^105 _ 0 


(79) 
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equations (38)^, (39)^, (40 ) k = 10> and (41) K = 1 „ give 

= J_ „3 10 - 5 ( Q! 7 Q! s + Q^a + «a«9)«in + 3(^7 + <*a + o?aMn - 2oj? n 

7106 60 10 a 6 (° ! 7 - a 6) (“8 - “a) (<*9 - “e) ~ 

(80) 

and correspondingly y 107 , 7 io 8> andy 109 . 

The first equations (42) and (43) yield 

y H3 = » yiu = - 774 - t 2- 784 • ( 81 ) 

c n c n c n 

From equations (38 ) k = 11> (39) (C _ 11 * ( 40 ^ k = ii’ and ( 41 ) K = n we can 
now obtain 


7117 = 


N(y 117 ) 

a 7 (a 8 - a 7 ) (a 9 - a 7 ) («io - « 7 ) 


(82) 


with 

N(y m ) = ^ a?i [10a 8 a 9 o; 10 - 5(a 8 a! 9 + a^ato + a 9 a 10 )« 11 + 3(a 8 + a 8 + aioMt - 2a?il 

- ril3 a 3( a 8 -«3)(° ! 9 -«3)( a 10 -“3)-ril4 Q! 4( Q! 8 - “ 4 ) ( a 9 -«4)( a 10 -“ 4 ) 


and corresponding formulas for y 118 , y 119 , andy 1110 . 

From the second equations (42) and (43) we find 

7 123 = “ T (c 7 773 + ^liy 113) > 7 124 = ~T (^7774 + ^8784 + ^Il7ii4) (83) 
u 12 M2 

For the computation of y 127 , y I28 , yi 29 , 7 i 2 io* and y 1211 we make use 
of (3S) K = i 2 , (39 ) k = 12 , (40) k _ 12 > (41 ) k = 12 and of 

C7Q75 + c 8 Q 85 + c 9 Q 95 + c 10 Q 105 + CjjQug + Ci 2 Q 125 = 
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The last equation is equation (IX, 20) ' of Table 1. If this last equation 
is satisfied, one of the error coefficients in x becomes zero, thereby- 
reducing the error in x . 

The above five equations yield 


y i27 


Q! 




O' 


NQymj 

7 )(«9 -« 7 )(“10 -“t)(«ii -«t) 


(84) 


with 

l 27 ) = 7 ^ [70q' 8 Q' 9 o' 10 o; 11 - 35(a 8 Q! 9 Q' 10 + o; 8 oi^a n + a z a io a n + O! 9 q' 10 o; 11 ) 
+ 21(q! 8 0' 9 + Q'gG! 10 + a % a n + 0' 9 0' IO + a 9 an + Q!ioaii) 

- 14(q! 8 + Q! 9 + a 10 + ail) + 1Q 1 
_ y 123 a 3 (a 8 - a 3 ) (a 9 -a 3 )(a 10 -a 3 )(a 11 - a 3 ) 

- Ti24 0! 4( a 8 - a i)( a 9 ~ a 4 )(ai 0 - a 4 )( a n - a \) 

~ 12 + ^ 3 ^ 6 " “ ^ 7< ^ 75 " ° 8< ^ 85 " ^ 9< ^ 95 “ ®ioQi 05 - C11Q115 ^ 


and corresponding formulas for 7 i 28 , yi 2 9 , Ti 2 io* an d Yi 2 n* 

The coefficients B and y can not be found from the equations of 
kO kO 

condition since they do not enter these equations. 

However, comparing (3) with (28), we find immediately 


B =01 - (B + B / 3 ) 

p k0 k V k1 h k2 m k,k-1' 


(k = 1,2,3 ...» 13) , (85) 


with the parentheses in ( 85) being omitted for k= 1. 



18. 


In No. 15 and No. 16 we have expressed the coefficients j3 and y 

kX kX 

by the coefficients a . As explained in No. 15, there are relations 
between the a These relations leave us with only four independent 
a : a 4 , a Jt a 9 , and o; 10 . Except for trivial restrictions, these four 
a can be chosen arbitrarily and lead to a Runge-Kutta-Nystrom formula 

K 

RKN-G-7(8) -13 based on thirteen evaluations per step of the differential 
equation (2) . 

19. There remains the problem of how to select the four independent coeffi- 
cients a 4 , a 7 , a 9 , ando! 10 . Naturally, one would like to have a Runge- 
Kutta-Nystrom formula with small truncation errors for whatever the 
problem (2) might be. Unfortunately, the truncation errors also depend 
on the problem (2) . A formula that might be very efficient for a certain 
problem can prove to be relatively poor for another problem. 

Therefore, the only reasonable way to select the independent parameters 
a 4 , a 7 , a 9 , q; 10 seems to be to find in Table 1 the error coefficients for 
the eighth-order terms in x and x . However, one has to keep in mind 
that these error coefficients have to be multiplied with certain expressions 
in the partial derivatives of ( 2) , summed up, and multiplied with h 8 to 
represent an approximation of the local truncation error in x or x . 

It can be assumed, however, that the local truncation error becomes 
small if the error coefficients are sufficiently small since the local 
truncation error obviously tends to zero if all error coefficients go to 
zero. 

From Table 1 we find the following 10 error coefficients in x which are 
different from each other: 

T rp rp rp rp rp rp rp rp rp 

16 » i 38> -*• 39 » 1 49> -*-67 » 1 78> i 81’ 1 82» 87 » 1 95 » 

the suffix indicating the number of the eighth-order equation of condition 
in Table 1. For instance 


c P — 

k k5 336 


12 

Tie = t 

K= 7 
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There are more error coefficients for x in Table 1. These additional 
error coefficients, however, differ by a constant factor only from those 
listed above. 

From the eighth-order equations of condition for x in Table 1 we find 
the following 28 error coefficients. 3 


T t , T 

16 » T 26 * T 48 ’ T 49 ’ 

T 59 > 

T 76 ’ 

T 89 > 

T 93 > T 

143 

Tl 59 > 

Tl75> 

^182 » 

Tl 89 > 

Tl 92 > 

Tl 94 > 

^213 * 

^220’ 

• • 

T 222 ’ T 223 

T 249 * 

^ 255 ’ 

^256 > 

^258 * 

T2G0’ 

^262 > 

^263’ 

^266 



It is, for instance 



K=7 


• 7 

c cr 

K K 


1 

8 


Again, there are more error coefficients for x in Table 1 which differ 
by a constant factor from those listed above. 

Since there are considerably more error terms contributing to the trunca- 
tion error in x than there are for the truncation error in x , the error 
control should be based on the truncation error in x . However, since 
this seems to be impossible without an unreasonable increase in the com- 
putational effort, we have to resort to an error control in x and try to 
keep the errors in x as small as possible. Since these errors in x 
propagate directly through the differential equation (2) , their influence 
is likely to be more serious than in the case of the differential equation 
( 1 ). 

20. We computed the above listed error coefficients in x and x for a large 
variety of combinations of the parameters a 4 , a 7 , a 9 , a 10 and finally 
decided on a combination for which the error coefficients in x as well as 
the ratio of the error coefficients in x to the error coefficients in x 
were reasonably small. 


3. Because of the choice of y 127 , . . , y 12 n in No. 16 the error coefficient 
^20 is zero. 
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We selected for our seventh-order Runge-Kutta-Nystrom formula RKN- 
G- 7(8) -13 the following combination 


1 3 1 3 , . 

«4 = 2 ’ “ 7 = 4 ’ = 8 ’ “ 10 = 8 ' ( 86 ) 

The coefficients for this formula are listed in Table 3. Since the last 
evaluation (28) is supposed to be taken over as first evaluation for the 
next step, the coefficients /3 130 » j 3 131 , , /3 1312 and y 130 , y 131 , .... , 

y 1312 in Table 3 have to be equal to the weight factors c 0 , c lf ... , c 12 
and c 0 , Cj, ... , c 12 of (29). 


The computation of the coefficients for our seventh-order formula and 
also for the sixth-order formula, listed later, was performed in 
40-digit arithmetic. 

Table 2 shows the pattern of our seventh-order formula. All coefficients 
different from 0 and 1 are marked by an asterisk. 


SECTION III. SIXTH-ORDER FORMULA RKN-G-617) 


21. The derivation of a sixth-order formula is similar to the derivation of 
the seventh-order formula in Section II. 

We base the sixth-order formula on ten evaluations with an eleventh 
evaluation which is taken over as first evaluation for the next step. 

Similar to Section II we make the following assumptions for the coefficients 
of our sixth-order formula: 


$1 = 

Cj = 0, 

ft - 

• • • » C4 — 

c 4 =0. c 5 = 

C5 , ... , 

Sb 

C 8 9 

c 9 = 

r . A _ 

C10 — C9 


Cl = 

0, . . . 

, c 4 = 0 ; 

«9= a 10 = 

1 






(87) 

*31 

-*41- 

= *ioi 

= ° Y 31 

= Y 41 = 

. = 

T 101 : 

= 0 




*52 

= *62 = 

* ‘ ' = *102 

= 0 *10 

,x = \ (A = 

= 0, 

1,2, . 

> . • 9 

9) 


(88) 




Y 10 

,x = 0 x (A = 

= 0, 

1,2, • 

1 • • y 

9) 




27 


and 


k— 1 


z wH (k=2 ' 3 - ••• • io) 

A— 1 


(89) 


K—l 

z 

\=1 


Vx-?**! <K = 2 ’ 3 ' ••• ’ I0) 


(90) 


K -1 

z 

X=1 


Vl = i“i (K = 5 ’ 6 10) 


(91) 


c 5^53 + c 6^63 + c 7^73 + c 8^83 + c 9^103 ® 

^ 5^53 + ^ 6^63 + ^ 7^73 + ^ 8^83 + ^ 9^93 = ® 


(92) 


c 5 / 3 54 + c e/ 3 64 + C 7^74 + c 8^84 + c 9^104 " 0 
® 5^54 + ® 6^64 + ' C 7^U + ^ 8^84 + ^ 9^94 = ® 


(93) 


and 


K — l 


S y KX 0l X" 6 a l (k " 2 ’ 3, "• ’ 10) 


A=1 


(94 


K-l 


2 vrs“l <K = 5 - 6 ' ••• • 10) 

A- 1 


(95) 


K — l 

z 

\=i 


v! = i”i (K ' 5 ' 6 10) 


(96) 
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c 5 t 5 2 + c eY62 + c 7'y 7 2 + c 8 y 82 + c 9 y 102 = 0 

c 5 y 5 2 + c 6 y 6 2 + c 7 y 72 + c 8 y 82 + c 9 y 92 = o 


( 97 ) 


22. The assumptions of No. 21 reduce the seventh- and lower-order equations 
of Table 1 to 

(11,1), (III, 1) , (IV, 1), (V,l), (VI, 1), (VII, 1), (VII, 10) 

The first six of these equations yield, in the same way as in Section II, 
the weight factors c 0 , c 5 , c 8 , c 7 , c 8 , c 9 and c 0 , c 5 , c 6 , c 7 , c 8 , c 9 . 
Equation (VII, 10) has to be solved together with the assumptions of No. 
21. The- resulting values for 0 21 , 0 32 , 0 42 , 0 43 , 0 53 , 0 54 , P 63 , 0 6 4> ^65 
and the restrictions for a lf a 2 , a 3 are the same as in Section II.. They 
are given by equations (46) through (53) . 

The remaining coefficients /3 , etc. are different from those of 

Section II. 

Putting 


B = 0 
73 


(98) 


we obtain from (89) k=? , (90) k=7 » (91) k _ 7 

J_ 2 . - 4(a B + a R )a 7 + 3a? 

74 12 7 a 4 (a 5 - a 4 )(a 6 - a 4 ) 


(99) 


and corresponding formulas for /3 75 and /3 76 . 
The first equations (92) and (93) yield 


^83 _ ( c 5^53 + c 6^63) 

^8 


^84 _ “ “ ( c 5^54 + c 6^64 + c 7^ 7 4) 

c 8 


( 100 ) 
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Equations (89) , (90) , (9l) , and(VII,10) can then be con- 

K — o K—o K — o 

sidered as four linear equations for the three coefficients j3 85 , j 3 86 , /3 g7 . 
They lead to a restriction for a s : 


a 8 


Nfa) 

D(« 8 ) 


( 101 ) 


with 

N(a 8 ) = 30(1 -a 5 )(l -a 6 )(l -a 7 ) [35a 5 o: 6 a 7 - 14(a 5 a 6 + cv 5 a 7 + a 6 a 7 ) + 7(o; 5 + a 6 + a T ) - 4] 
+ [30a 5 Q! 6 Q: 7 - 20(a 5 o- 6 + a- 5 or 7 + a 3 a 7 ) + 15(a 5 + a 6 + a 7 ) - 12] • 

[35o- 6 a 6 a 7 - 2l(a 5 a: 6 + a- 5 a 7 + a- 6 a- 7 ) + 14(a 5 + o- 8 + a- 7 ) - 10] 

D(a 8 ) = 30(1 -ar 5 )(l -a e )(l -a 7 ) - 14 ( 0 - 50-8 + a 5 o- 7 + a- G a 7 ) + 7(ce 5 + a 6 + a 7 ) . 4] 

+ [SOajQ^o-j - 20(a- 5 a- 6 + a s a 7 + a e a 7 ) + 15(a 5 + o: 6 + a r ) - 12] • 

[70a 5 a 8 a- 7 - 35 (o- 5 q- 6 + o- 6 o- 7 + a 6 o- 7 ) + 21 (q- 5 + a- 6 + a 7 ) - 14] 


and to 


086 


jj a|[6o- 6 a 7 - 4(a 6 + a 7 )a s + 3ajj] -0 83 a 3 ( - a 3 )(a 7 - a 3 ) -0 84 a 4 ( Q! 6 - a t )(a 7 - o- 4 ) 


( 102 ) 


and corresponding formulas for /3 86 and j 3 87 . 
The second equations (92) and (93) give 


^93 T~ (C 5 P 53 + ^63 + ^ 8^83 ) 

c 9 

P 94 = “ JT (^ 5^54 + ^ 6^64 + c lPu + ^ 84 ) 
c 9 


(103) 


30 



095 ~ 


23. 


Equations (89) , (90) , (9l) , and(VII,10)‘ then represent 

k — y k — y k — y 

four equations for the four coefficients /3 95 , j3 97 , /3 98 . Their solution 

is 

1 lOa^aj - 5(0^ + a 5 Q!7 + <**<* 7 ) + 3(« 5 + «« + « T )-2 , 1ft4 v 
^ 98 60 c 9 a 8 (a 7 l"a 8 ) (a 6 -a 8 )(o! 5 - a 8 ) ' ’ 


[6a 6 Q! 7 - 4(a e + a 7 ) + 3] -^ w of 3 (o( 6 


-«3)( a 7 -«3) -094“4(«6 " “ 4 ) («7 - “ 4 ) “ ^98«s( Q! e - “fl) ( a 7 - « B ) 

a 5 (a 6 - a 5 )(«7 - ^ 5 ) 


(105) 


and expressions /3 9g , /3 97 that correspond to /? 95 . 

We still have to determine the coefficients y . The coefficients y 21 , 

y 32 , T 42 > Y 43 » Y 5 2 > Y 53 . Y 54 » Ys 2 » y 6 3 > Yg 4 > Yes are the same as the corre- 
sponding coefficients in Section II and are given by equations (68) through 
(73). 

Setting 


Y72 - Y 7 3 _ 6 , 

we find from (94) (95)^, (96) k=7 


Y 74 


J_ a z 10 Ml: 5 ( a s + ggj 

60 7 a 4 (a! 5 - « 4 )(a 6 - of 4 ) 


and corresponding expressions for y 75 and y 76 . 
The first equation (97) yields 


Y82 — “ Y 52 

c 8 

Setting 

Ys 3 = Y84 = 0 


(106) 


(107) 


(108) 


(109) 
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equations (94) K=g , (95) K=g , (96) K=g give: 

1 10a 8 a 7 - 5(a e + a 7 )a 8 + 3oi\) -y i2 o> 2 {a 8 - a 2 )(a 7 - oi 2 ) 
* ~ a 5 ( a 6 -“5) (“7 -“5) 

and corresponding expressions for y 86 and y 87 . 

From the second equation (97) , we find 


(110) 



(^5^52 + ^8^82) 


With 


( 111 ) 


T93 _ T94 - 795 _ 0 


( 112 ) 


we find from (94) k=9> (95) K=9 * ( % ) K=9 


[ 10o; 7 a 8 - 5(a 7 + a 8 ) +3] - y 92 a 2 (a 7 - a 2 ) (a 8 - a 2 ) 


y% = 


« 6 (« 7 -«e)(" 8 - "e) 


(113) 


and corresponding expressions for y 97 and y 98 . 

This concludes the computation of the coefficients |3 and y , since 

K A K A 

(3^ and y are again obtained from (85) with k = 1,2,3, ... , 10. 

24. The expressions for |3 and y^ in No. 22 and No. 23 contain four 
parameters a 4 , a 5 , a 6 , and a 7 which we can choose arbitrarily. 

As in Section II, we choose these parameters so that the error coefficients 
in x and the ratio of the error coefficients in x and in x become 
reasonably small. 

In the case of a sixth-order Runge-Kutta-Nystrom formula there are five 
different error coefficients in x: 


T rp rp np rp 

10* a 23> i 24» i 28> 1 34 
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and eighteen different error coefficients in x: 


T rp rp rp rp rp rp rp rp rp 

1» J- 10» M2> i 16» i 31> i 32» 38 » 1 39> A 47* 1 49 

• ■ • • • « • • 

T rp rp rp rp rp rp rp 

81* i 82» i 86> -*-87* -*-89* i 91* i 92» *-95 


In our final sixth-order formula we chose for the free parameters 


OL, 


5 

8 ’ 




(114) 


The coefficients of the sixth-order formula with the parameter values 
(114) are listed in Table 5. 

The pattern of our sixth-order formula is shown in Table 4. 


SECTION IV. FIFTH-ORDER FORMULA RKN-G-5(6) 


25. Although it is possible to construct fifth-order formulas based on seven 
evaluations per step, we prefer to use eight evaluations per step, since 
we then obtain formulas with smaller local truncation error terms. In 
spite of the one additional evaluation per step such formulas proved to be 
more economical. They allow a larger stepsize because of their smaller 
truncation errors. 


We make the following assumptions for the coefficients of our fifth-order 
formula: 


II 

<C? 

o 

II 

T-« 

o 

II 

rH 

<o 

c 2 0 , 

A 

o 

CO 

II 

O 

5 

A 

A 

A 

A 

= 0 , 

A 

O 

II 

o 

45*. 

c 5 “ c 5 > 

c 6 - c 6 

, c 7 

Cs - c 7 

Ci = 0 , 

c 2 = 0 , 

c 3 = 0 , 

a 7 = 

8 

CO 

II 

I- 1 



(115) 
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031 " 041 _ 051 " 061 “ 071 ” 081 ” 0 
741 = 751 = 761 = 771 = 781 = 0 
0 8A = c x (*= 0,1,2, ... , 7) 

7g^ = (X= 0,1,2, ... ,7) 


and 


(116) 


K-l 

E 

A=1 


B a = — 
h kA X 2 


cr 


(« = 2,3, 


8 ) 


(117) 


K-l 

Z 

\=i 


0 

kX X 


2 _ 


= — 01 
3 K 


( k — 2,3, 


8 ) 


c 4042 + c 5052 + c 60 62 + c 7082 ~ ® 
^4042 + ^5052 + ^6062 + ^70 72 = ^ 


and 


K-l 

E 


\=1 


7 

kX X 



(k = 2,3, 


8 ) 


(118) 


(119) 


( 120 ) 


K-l 

E 

\=1 


y a 2 ~ 
y kX X 


= J2 a -' (K= 2 ’ 3 ’ 


8 ) 


(121) 


^4742 + ^5752 + ^6762 + ^7772 - ® • ( 122 ) 

26. The assumptions of No. 25 reduce the sixth- and lower-order equations of 
condition of Table 1 to 

(11,1), (111,1), (IV, 1 ), ( V, l) , (VI, 1 ), (VI, 6) 
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The first five of these equations yield, as in Section II and III, the weight 
factors c 0 , c 4 , c 5 , c 6 , c 7 and c 0 , c 4 , c 5 , c 6 , c 7 . Equation ( VI, 6) has to 
be solved together with the assumptions of No. 25. The resulting values 
for /3 21 , 0 32 and the restrictions for a 2 are the same as in Sections 
II and III and are given by equations (46) through (49) . 

Equations (119) can be satisfied by 

^42 = @52 = ^62 = ^72 = @52 = 0 • ( 123) 

From(ll7) , and ( 118) „ we then find 

' 7 k=4 k=4 

a 3 = | a 4 (124) 

and 

^43 = ^ a 4 • (125) 


Equations ( 117) and ( 118) yield 
K — 5 K — D 


@53 


1 9 3 Q! a — c 

— ai — -A \ 

6 5 a 3 (a 4 -a 3 ) 


£54 


1 o 3a 3 - 20! = 

— G!c — 7-“ \ 

6 5 G! 4 (0!3 - a 4 


Equations (ll7)^_g, (118)^ and (VI,6) result in 


^63 


J_ a 2 - 4(oii + a B )«B + 3QigRi 

12 6 4 - a 3 )(o; 5 - a 3 ) 


(126) 


(127) 


and corresponding expressions for /3 64 , j3 65 . The abbreviation R 4 in 
( 127) stands for 


Ri 



C 4 / 3 43«3 - C 5 (@53 a 3 + / 3 54«D " J ° 7 
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We compute the coefficients j3 73 , j3 74> j3 75 , /? 76 from(ll7) , (118) , 

(VI,6) ' and from the additional condition T l0 = 0 (see Table l) . The 
last condition helps to reduce the local truncation error in x. The above 
equations lead to 


= J_ 30 OjOja^ - 20(o; i Q; s + a 4 o? fi + 01 ^) + 15^ + a, + a fi ) • R ? - 12R 3 
73 6 0 ff 3 («4 - «3)(“5 - «3)(° ; 6 - “a) 

(128) 

and corresponding expressions for /3 74 , /3 75 , j3 76 . Here we have used the 
abbreviations 

R 2 = ^ * c 4^43 o; 3 - ^(/^af + / 3 54 0? 1) 

R 3 = ^ - c 4 /3 43 a| - c 5 (/3 53 a| + /3 54 a 4 ) 


- c 6 (/3 63 a| + /3 64 0!^ +/3 65 Q! 3 5 ) 

“ ^(^63® 3 + ^G4 a 4 + ^65 a 5) 


27. We now compute the coefficients y . The coefficient y 21 is again given 

kA 

by (68). The coefficients y 3l and y 32 are obtained from equations 
(120) K=3 aad(121) K=3 : 


Tsi “ 77 “s 


2Q! ? - ou 


12 3 o' 1 (a 2 - a 4 ) * 


T32 


= J_ a 3 2a t -a 3 


12 3 a 2 (a 4 - a 2 ) 


(129) 


With y 41 = 0 we find from (120)^_^ and (l2l)^_^t 


1 o 2a 3 - a 4 

742 = 12 *'* ^ 


12 " 4 a 3 (a 2 - a 3 ) 


3 2 ^2 - a 4 


(130) 


Setting 


^62 _ '>'72 " 0 


(131) 
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we find from ( 122) : 


c 4 

752 = - 7 T ?42 

Equations (120) and ( 121) then yield 

K — O K — 0 



With 


(132) 


(133) 


! 

763 = r?3 = 0 


(134) 


we find from ( 120) and(l2l) . : 

' ' k=6 k=6 

1 , 2o? C - a R 1 

= 12 “« ' T6i = 12 

and from (120) , (121) _ and from T 8 = 0 (see Table l) : 

J_ 10a s «fi - 5(0=5 + «r) + 3R 4 
774 60 a4(a 5 - a 4 ) (a 6 - “ 4 ) 



(135) 


(136) 


and two corresponding formulas for y 75 and y 76 . Here we have used the 
abbreviation 

r 4 = - c 4 y 43 a| - c 5 (y 53 o!l + y 54 o;|) - 6 6 (y G3 G!| + y u a\ + y 6 5 «i)J . (137) 

This concludes the computation of the coefficients /3 and y ^ since 

/3 and y are again computed from (85) with k= 1,2,3, ... , 8. 

K 0 K 0 
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28. The expressions for )3 and y ^ in No. 26 and No. 27 contain three 

free parameters: a 4 , a 5 , and a 6 . By a proper choice of these param- 
eters, we might try to obtain formulas with reasonable small error 
coefficients in x as well as in x . 

In the case of our fifth-order formula we have only one error coefficient 
in x: 


t 6 


and four error terms in x: 


Ti> Tg, T 24 , T 34 

29. It is interesting to notice that by adding one more condition ( VII, l) ' to 

(II, l) ’ , (lH,l)’, (IV, l) ’ , (V,l)’, and ( VI, l) * , we obtain a restrictive 
condition for the a 1 s: 

5a,ac - 3 (q! 4 + Q!r) + 2 
01 6 ~ 10 a 4 a 5 ' - 5(a 4 + a 5 ) +3 

that makes all our error coefficients T 6 , Tj, T g , T 24 , T 34 zero. Natu- 
rally, such a choice of a 6 is not suitable for our fifth-order formula, 
since our stepsize control would break down in this case. 

However, by choosing a 6 close to the value of (138) , we can obtain 
sufficiently small error coefficients in x and x that lead to efficient 
fifth-order formulas. 

Such a formula is obtained for 



a 4 = 


_9_ 

10 ’ 




(139) 


For the values cn 4 and a 5 of (139) the condition (138) would result in 
a 6 = 17/60 (« 0. 2833) , which is reasonably close to a 6 = 2/7 (« 0. 2857) . 
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The coefficients for our fifth-order formula based on the values (139) for 
the free parameters c^, as, a 6 are listed in Table 7. Since these 
coefficients have relatively simple values, we have listed them in frac- 
tion form. 

The pattern of our fifth-order formula is shown in Table 6. 


SECTION V. APPLICATION TO TWO NUMERICAL EXAMPLES 


30. In this section we apply the Runge-Kutta-Ny strom formulas of this report 
and the Runge-Kutta formulas of [3] to two problems; one of them is 
linear in the first derivatives and the other one is nonlinear. As the 
linear problem we choose an orbit of the restricted problem of three 
bodies. The same orbit has already been integrated in an earlier paper 
of ours [5] , using a power series expansion technique. The orbit is 
pictured in Figure 1. 

Table 8 shows the differential equations and the initial conditions for the 
problem. Since the problem has no solution in closed form, we integrated 
the problem by the above mentioned power series expansion technique 
using thirty decimal digits. 


Truncating the series after 12th-order, 16th-order, or 20th-order terms, 
the results (for t= 6)for x, y, x , y agreed to about twenty decimal 
places. Rounding these results to 16 decimal places, we found 


t = 6 


0. 1167 

0361 

7342 

5520 

0. 1966 

6280 

9565 

9560 

0.3341 

2960 

6037 

2482 

0.9745 

3805 

7977 

8027 


(140) 


We substituted these values for the solution of our problem. The errors 
Ax, Ay, Ax, Ay in Table 8 are the deviations of the solution obtained 
by our Runge-Kutta -Nystrom or Runge-Kutta formulas from the above 
values ( 140) . 
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31. Table 9 shows the differential equations and the initial condition for a 
problem that is nonlinear in the first derivatives. Since this problem 
has a solution in closed form, the errors of our numerical solutions 
could easily be established. 

32. All calculations in Tables 8 and 9 were executed on an IBM -7094 computer 
in double precision (16 decimal places) . The computer was equipped 
with an electronic clock to measure the execution time for the various 
formulas. 

The stepsize control for our Runge-Kutta-Nystrom and our Runge-Kutta 
formulas is described in No. 26 of our earlier report [l] . 

33. Tables 8 and 9 show the results of the various formulas applied to 
Problem I and II. Comparing the results of our new Runge-Kutta- 
Nystrom formulas with those of our Runge-Kutta formulas of [3] , we 
notice that, in these examples, we save from 25 percent to 60 percent 
of the execution time by using the new formulas. In most cases, our 
new formulas are also slightly more accurate. When relaxing the toler- 
ance from 0.1 • 10 -16 to 0.1 • 10 -15 , the relative savings in execution 
time for our new formulas do not change very much. 


Computation Laboratory 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama 35812, May 31, 1974 
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TABLE 1. (Continued) 
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TABLE 1. (Continued) 
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TABLE 1. (Continued) 



/k-1 

E c a* E * a P 

^ K K \ “ K\ X \1 
3 \ 2 * 


/k-1 

y c a 2 y p a 
3 * * \ 2 * 


2 c r a *r ( 2 ^ P > 


« k k\ U X2 

3 V 2 


/k-1 > 

2 7« ( 5 V*A P A1 

3 \ 2 } 


/k-1 

2 °k“k ( 2 \a Q A1 

3 \ 2 


= ? V«(? T ' 


| V« (2 *«a“x P A1 


= y, c fl y )3 ,a^Q. 

“ K K \ „ K\ X , 


3 \ 2 


/k-1 

S YM £ /3 KX a \ P A2 


/K-1 > 

3^ Ck ° \ 2^ ^ kxQ ^ 2 


E Z ^«-A P A 


o K K \ V r d A3 

3 \ 2 




y c f y 1 y 
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TABLE 1. (Continued) 
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X 


2 


c , c 


vin 


/k-1 ' 

2 o K ( S ^X P U 

3 \ 2 


Z c Ai(pKx p xi 


/k-1 

Z C k P k2 ( Z P ** P M 
3 \ 2 


/k-1 

2 C K P Kl ( ^ y KX P Xl 

3 \ 2 


y c p ( y t$ a p 

^ K Kll “ KA A A 1 
3 \ 2 , 


/k-1 

? Wl 2 

3 \ 2 


/k-1 

2 c « p JpA 

3 \ 2 


K-1 / A-l > 

£ c k“k 1^(1 Wl 


K-1 /A-l 

Z c K a K Z Y K xy Z l 


1680 Z c K a K Z ^“xl Z 


K-1 /A-l 


4 3 


K-1 /A-l 

Z C K a K Z ^A^ Z Y XM P /il 


K-1 /A-l 1 

Zca Z^i( Z ^ P 1 

^ k k 3 K A \ 2^ Am M Ml 


K-1 /A-l 

J C k“k pKX^W 

K-1 /A-l N 

4 ° K “ K 3 M I 


J C K | y K\“x( | VV 


65 


1 

10 080 


1 

1260 
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TABLE 1. (Continued) 



i 
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TABLE 1. (Continued) 







X 


/K-l 

^ C K Q! K P Kl ( 2 Vx P Xl 

3 \ 2 


/K-l 

^ WkI ( ^ /3 kX P X2 
3 \ 2 , 


S °K ( | \x a i p xi 

3 \ 2 


I C \t 7 ^ 1 


? C| < ( ? 7 ^ Px2 


/k-1 

^ C k ( ^ Y K\“x Q X 
3 \ 2 


5 m 5 Ykx “ xPx3 




£ C K f ? \aSu 


? 


/ K-l 

X C K ( 5 7 kX P X4 

3 \ 2 


/k-1 

^ 1 ^ 7 kA P X2 P X1 

3 \ 2 y 


5*- (M (?'•'» 


? 4‘t ** a * p ™ 

3 \ 2 


VID 









TABLE 1. (Continued) 
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95 


96 


97 


98 


99 


100 


101 


102 


103 


104 


105 


106 


107 


108 


1 

4536 


I* C * I? 


1 

504 


VIII 


1 

4032 


10 368 


1 

3456 


£ C k ( ^kX P X1 
3 \ 2 


?'■ 


K-l \ / K— 1 

2 ^kX P Xi) ( / ^kxS 1 


K-l \ /K-l 

2 ^KX P Xlj ( ^ ^kx“x P X1 


— = 7 

5184 u 


K-l \ A-1 

? ( ? 


B P 
kX X2 


5184 ? C k^k2 ( ^ ^kX P X1 


10 368 


/ K-l 

^ C K Q Kl( ^ Y kX P X1 


1 

448 


1 

1152 


1 

384 


1 

576 


1 

576 


1 

1152 


1 

3456 


10 368 


K-l 


£ c k Q ki( £ ^ K \ a \ p ; 
3 \ 2 

S C K Q Kif ? ^KX Q X1 


XI 


1 

384 


1 

1152 


1 

5184 


| C k Q k1 ( | ^ P X2 


1 

576 


1 

1728 


^ °k P k3 ( ? ^kX P X1 


1 

192 


1 

5184 


/K-l 

^ C k P k2 ( ^ Y kX P X1 
3 \ 2 


1728 " ^ °K P K2 ^ | 0 kx“\ P M^ 


1 

5184 


1 

2592 


A-i 

^ C k P k 2 y ^ ^kX Q X1 

/ K-l 

S C k P k2\^ ^ /3 kX P X2 


1 

576 


1 

192 


1 

576 


109 


1 

288 
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! 

f 
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X 


X 


K— 1 A-l /m-1 > 

5^ | ^ ^ KX 3^ 2^ y ^ Vvl i 

Ik- 1 pul h-1 

JV*K \Z*<K 


f K~1 f A-l /fl-1 \1 I 

8 e ' t “ K | 4 ^ \ M I ^^■'7 ) 


K— 1 A-l A*-l > 

Z °k“k 2 P K\ 2 ^(-? V P I>2 

5 ^ 4 |_ 3 \ 2 y 

f K-l A-l /ju-1 

2 °k 2 y kx"\ 2 V V p »i 


K-l A-l Ai-1 

Z c j Z Va Ev Sl»„ 


5 K 4 kA 3 X ^\ 2 ^ Vl 


K-l A-l Ai-1 

E c J Z v KX E p x « ( E <** 


5 K 4 kX 3 XM M \ 2 Vl 


K-l A-l Ai-1 
^ C k| ^ 7 kA ^ ^A/i ( ^ y /ii/ P i/l 


K-l A-l /jU-1 

^ °k| £ 7 kA ^ *Vy ^ 

( K-l I" A-l /m- 1 \ 

2 C K | | 7 ka | | V Q >*J 

f K-i r A-i /\i~\ \ 

S C K | ^ y ka ^ ^ 

2 J Jf.,.; [ zV/'z * p . 


^ K ^ kA A 3 Am l 2 vl 


K-l A-l /m-1 

^ C K ^ ^K\ a \ ^ ^ \ ^ vl 


K-l A-l At- 1 

E c K | E V*x E VVy | V p *i 


K-l A-l Ax-1 

^ C K ^ P KA a A | 2 


vm 


229 


1 

60 480 


1 

6720 
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TABLE 2. PATTERN FOR RKN-,G-7( 8) -13 


0123456789 10 11 12 0 123456789 10 11 12 


0 0 0 


* 0 * 


i * o * 


* o * * 


1 * 00 * 


* 00 ** 


^ o ^ ^ ^ 


*00*** 


* 0 0 * * * 


* * 0 0 0 0 ** 


* 0 0 * * * * 


8 ** 0 0 0 0 *** 


* 000 **** 


9 * * 0000 **** 


* 0000 **** 


10 ** 00000 **** 


* 00000 **** 


* 0000 ***** * 


* 00 ** 00 *** * 


* 0000 ***** * * 


* 00 ** 00 *** * * 


* 000000 *** * * * 1 * 000000 *** * * * 


n 
















TABLE 3. COEFFICIENTS FOR RKN-G-7( 8) -13 


Of 1 


0. 7347 

0804 

8410 

643 8 

3606 

1035 

6687 

4183 • 10" 1 

O' 2 

= 

0.1102 

0620 

7261 

5965 

7540 

9155 

3503 

1127 

«3 

= 

0.1653 

0931 

0892 

3948 

6311 

3733 

0254 

6691 

O' 4 

= 

0.5 








Of 5 

= 

0.2662 

8826 

9291 

2616 

4263 

5203 

6943 

9706 

CVg 

= 

0.6337 

1173 

0708 

7383 

5736 

4796 

3056 

0294 

O' 7 

=r 

0. 75 








0'S 

= 

0. 5625 








O' 9 

= 

0.125 








O' 10 

as 

0.375 








O' 11 

= 

0.9665 

2805 

1602 

3570 

0834 

4516 

4211 

9099 

O' 1 2 

= 

1 








O' 13 

= 

1 









010 

= 

0. 7347 

0804 

8410 

643 8 

3606 

1035 

6687 

4183 

• 

10 _1 

020 

= 

0.2755 

1551 

8153 

9914 

3852 

2888 

3757 

7819 

• 

10" 1 

0 21 

= 

0. 8265 

4655 

4461 

9743 

1556 

8665 

1273 

3456 

• 

10 -1 

0 30 

= 

0.4132 

7327 

7230 

9871 

5778 

4332 

5636 

6728 

• 

10 _1 

0 31 

= 

0 










0 32 

= 

0.1239 

8198 

3169 

2961 

4733 

5299 

7691 

0018 



040 

— 

0. 8967 

0558 

7637 

9578 

2658 

3783 

2538 

9875 



041 

= 

0 










042 

= 

-0.3458 

5915 

2683 

1433 

6799 

4110 

9332 

7799 

• 

10 +1 

0 43 

= 

0.3061 

8859 

3919 

3475 

8533 

5732 

6078 

8811 

• 

10 +I 

050 

= 

0. 5705 

3369 

4239 

6532 

8229 

3636 

4664 

4660 

• 

10 _1 

051 

= 

0 










052 

~ 

0 










053 

= 

0.2066 

4670 

6954 

9824 

5793 

1750 

7585 

7561 



054 

= 

0.2588 

1929 

1231 

3856 

4740 

8928 

9176 

7848 

• 

CM 

o 

iH 

060 

= 

0.2213 

0953 

4027 

3273 

8534 

2980 

542 8 

5043 

• 

10 _1 

061 

= 

0 










062 

= 

0 










063 

= 

0.3666 

6901 

8421 

5938 

0713 

9353 

1564 

8106 



064 

= 

0.3207 

5560 

5327 

6707 

8702 

4980 

3819 

9180 



065 

= 

-0. 7584 

3846 

4432 

5897 

5333 

8352 

8715 

4967 

• 

10 -1 

0 70 

= 

0. 8333 

3333 

3333 

3333 

3333 

3333 

3333 

3333 

• 

10 -1 


071 - 0 

07 2 ~ 0 

073 = 0 

74 = 0 





TABLE 3. (Continued) 


0 75 


0.3843 

6436 

9641 

3162 

1037 

9110 

0848 

8971 


076 

= 

0.2823 

0229 

7025 

3504 

5628 

7556 

5817 

7696 


0 80 

= 

0. 8349 

6093 

7500 

0000 

0000 

0000 

0000 

0000 • 

10 -1 

@81 

= 

0 









@82 

= 

0 









0 83 

= 

0 









P 84 

= 

0 









P 85 

= 

0.3830 

6747 

4793 

9740 

8845 

8700 

6356 

1136 


086 


0.1354 

8721 

2706 

0259 

1154 

1299 

3643 

8864 

10" 1 

087 

= 

-0.3955 

0781 

2500 

0000 

0000 

0000 

0000 

0000 • 

@20 

= 

0. 7342 

0353 

2235 

9396 

4334 

7050 

7544 

5816 • 

10" 1 

P 91 

= 

0 









@92 

= 

0 









@93 


0 









@94 

= 

0 









@95 

= 

0.9880 

8964 

9160 

2291 

6024 

2057 

1033 

64 20 • 

10 _1 

@98 

= 

0.2415 

3311 

3273 

2774 

9549 

84 2 8 

0345 

1955 

10 _1 

@97 

= 

-0.4870 

7561 

7283 

9506 

1728 

3950 

6172 

8395 • 

@98 

= 

-0.2400 

54 86 

9684 

4993 

1412 

8943 

7585 

7339 


@100 

= 

0. 8137 

8441 

1270 

6706 

4904 

0410 

5640 

4207 • 

10" 2 

@ 101 

= 

0 









@102 

= 

0 









@103 

= 

0 









@104 

— 

0 









@105 

= 

0 









@106 

= 

-0.3626 

6091 

1746 

4713 

4384 

0315 

3205 

8792 

10 _1 

@107 

= 

0.6972 

6880 

5971 

2792 

8317 

2726 

0984 

7243 • 

@108 


0.3779 

7780 

6207 

6339 

2161 

1543 

4150 

9711 


@109 

= 

0.2818 

1838 

0829 

0027 

8742 

1095 

1900 

0315 

+ 
o 
1 — / 

@110 

= 

-0.1404 

2538 

9224 

8283 

8913 

2800 

3122 

54 7 6 • 

@111 

= 

0 









@112 

= 

0 









Pm 

= 

0 










0114 = 

0 


0115 = 

-0.1355 

5559 

0116 = 

-0.1502 

1472 

0117 = 

0.1476 

7543 

0118 = 

-0.2170 

7681 

0119 ~ 

0.6614 

9759 


0294 

0495 

7528 

3041 

8248 

4805 

0961 

7213 

2841 

6794 

9686 

2336 

9651 

3368 

8432 

5773 

5026 

7655 

8681 

0392 


1334 

2361 * 

10 +2 

3096 

9968 • 

10 +1 

0684 

1588 • 

10 +1 

7360 

7995 • 

10 +1 

0283 

3030 • 

10 +1 


72 




TABLE 3. (Continued) 


Pino 

- 

0.1150 

7526 

1735 

6932 

1530 

6792 

2237 

6434 

• 

10 +2 

0120 


- 0. 5270 

8651 

8158 

0131 

5268 

1768 

8218 

7497 

• 

10 +1 

0121 

= 

0 










0122 

=: 

0 










0123 

= 

0 










0124 

= 

0 









10 +2 

0125 


- 0.4996 

5599 

5536 

5683 

3001 

0459 

2152 

9105 

• 

012 S 

= 

- 0. 5030 

2228 

9286 

5823 

1516 

1351 

2481 

2231 

• 

10 +1 

0127 


0.4454 

8269 

0452 

9876 

0506 

5182 

3862 

2704 

• 

10 +1 

0128 

= 

- 0. 8607 

1533 

1240 

3384 

1312 

4067 

4298 

9148 

• 

10 +1 

0129 


0.2384 

0410 

0463 

7228 

7590 

0786 

7645 

6468 

• 

10 +2 

01210 

= 

0.4171 

1581 

4660 

2838 

8124 

0696 

6716 

4840 

• 

10 +2 

01211 

= 

- 0.1329 

7747 

6424 

3799 

54 0 8 

2370 

9555 

8512 


10" 1 

0130 

= 

0.3509 

9303 

0565 

8188 

3152 

6601 

7368 

1744 

• 

0131 


0 










0132 

= 

0 










0133 

= 

0 










0134 

— 

0 










0135 

= 

0 










0 136 

= 

0 










0137 

= 

0.2522 

3475 

2766 

3160 

6400 

6388 

5341 

7712 



0138 

= 

0.1184 

0033 

3068 

7654 

9234 

1625 

1536 

4336 



0139 

= 

0.2025 

8133 

6112 

5092 

9893 

1878 

9987 

1888 



01310 

= 

0.2675 

7025 

2594 

2014 

0796 

3933 

2927 

2621 



01311 

= 

0.1658 

6384 

5106 

2987 

3791 

2680 

9815 

0965 



01312 

= 

- 0.4174 

9822 

7046 

7288 

4309 

1671 

3445 

6960 

« 

10" 1 

7io 

__ 

0 . 2698 

9795 

8199 

6884 

8329 

9949 

7050 

8715 

• 

10“ 2 

Y20 

= 

0.3036 

3520 

2974 

6495 

4371 

2443 

4182 

2304 

* 

10" 2 

Y21 

= 

0.3036 

3520 

2974 

6495 

4371 

2443 

4182 

2304 

• 

10" 2 

Y30 

= 

0.6831 

7920 

6692 

9614 

7335 

2997 

6910 

0184 

t 

10“ 2 

73 1 

= 

0 









10“ 2 

732 

= 

0.6831 

7920 

6692 

9614 

7335 

2997 

6910 

0184 

m 

740 

= 

- 0.1026 

3757 

7319 

7788 

8994 

3108 

7282 

4217 

• 

10 " 2 

741 

= 

0 










742 


0 










743 

= 

0. 1260 

2637 

5773 

1977 

8889 

9431 

0872 

8242 



750 

= 

0.9890 

9903 

8431 

0741 

7913 

3134 

9906 

4241 

• 

10 -2 

751 

= 

0 









10" 1 

752 

= 

0.2040 

1758 

7591 

1134 

9514 

1705 

1849 

8571 

• 

753 

= 

0 . 5026 

5147 

7133 

2870 

3261 

8251 

0473 

5338 

• 

10“ 2 


7 


73 




TABLE 3. (Continued) 


T 54 = 

0.1354 

5726 

6312 

7775 

5415 

7283 

6001 

4730 

1 

10 -3 

760 = 

0.3677 

2464 

6953 

1772 

1429 

7415 

7246 

2201 

• 

10" 1 

Yei = 

0 










762 = 

0 










763 = 

0. 8213 

2294 

7785 

2178 

5827 

7217 

4140 

7693 

• 

10 -1 

764 = 

0.3008 

7165 

4090 

9896 

3036 

8709 

1811 

9641 

• 

10 _1 

765 = 

0. 5180 

3353 

9359 

9379 

0519 

8241 

0553 

1789 

• 

10 -1 

770 = 

0.4123 

3049 

0882 

7287 

3123 

2210 

0402 

1091 

• 

10 -1 

771 = 

0 










ro 

II 

0 










773 = 

0.1133 

5100 

2930 

6181 

9105 

3287 

9807 

8376 



774 " 

0. 5672 

2148 

5922 

3766 

8841 

3016 

7743 

6715 

• 

10" 1 

775 ~ 

0. 5745 

6202 

0649 

54 52 

5469 

3769 

2473 

6474 

• 

10 _1 

776 = 

0.1248 

7597 

3239 

1674 

1512 

8124 

1302 

1961 

• 

10 _1 

780 = 

0.4214 

6301 

2695 

3125 

0000 

0000 

0000 

0000 

• 

10 -1 

781 = 

0 










782 = 

0 










783 = 

0 










784 = 

-0. 7808 

8073 

7304 

6875 

0000 

0000 

0000 

0000 

• 

10" 1 

785 = 

0.1410 

4682 

1029 

2877 

2004 

3970 

8553 

6135 



786 = 

0. 7460 

3813 

7363 

3727 

9956 

0291 

4463 

8648 

• 

10 _1 

787 = 

-0.2150 

5737 

3046 

8750 

0000 

0000 

0000 

0000 

• 

10 _I 

790 = 

0. 5524 

3877 

1719 

2501 

1431 

1842 

7069 

0444 

• 

10' 2 

791 = 

0 










792 = 

0 










793 = 

0 










794 = 

0 










795 = 

0.4591 

3375 

8935 

0515 

8838 

0181 

1180 

7029 

• 

10~ 2 

796 = 

0.1200 

9956 

9922 

6813 

9808 

9279 

5562 

3138 

• 

10" 1 

797 = 

-0.2436 

1818 

4156 

3786 

0082 

3045 

2674 

8971 

• 

10" 2 

798 = 

-0.1187 

7000 

4572 

4737 

0827 

6177 

4119 

7988 

• 

10 _1 

7ioo = 

0.1239 

6099 

0923 

0042 

8073 

8555 

8121 

1091 

• 

10 _1 

Y101 = 

0 










7102 = 

0 










7103 = 

0 










7104 = 

0 










7105 = 

0 










7106 = 

-0.2314 

8568 

8348 

8114 

9606 

8286 

3748 

4335 

• 

10 _1 

7 1 07 = 

0.4405 

7716 

3385 

9229 

4670 

5995 

3820 

0368 

• 

10" 2 

7108 = 

0.2416 

4236 

8703 

9608 

6181 

8918 

2892 

7171 

• 

10' 1 


74 



TABLE 3. (Concluded) 


Yl09 

= 

0. 5249 

4961 

2383 

2540 

5884 

0212 

7352 

6037 

• 

10" 1 

TllO 

= 

-0.1214 

8292 

3371 

7236 

6838 

6923 

7170 

6654 



Till 

= 

0 










7112 

= 

0 










7113 

= 

-0.1594 

8786 

8094 

6904 

7245 

6588 

6876 

3595 

• 

10 +1 

7114 

= 

0.7708 

9844 

4095 

9035 

4601 

1435 

8063 

0038 

• 

10' 1 

7115 

= 

0 










7116 

= 

0 










7117 

= 

0.9884 

4932 

1354 

4261 

8048 

6243 

2844 

1900 

• 

10"" 1 

7118 

= 

-0.1851 

7690 

1776 

5400 

9760 

1245 

5930 

3975 



7119 

= 

0.1666 

5727 

1178 

0734 

2381 

8671 

54 63 

0279 

• 

10 +1 

71110 

= 

0. 5261 

1925 

5036 

5255 

2568 

0405 

9902 

2802 



7120 

= 

-0.4947 

5846 

7641 

0233 

2689 

6037 

0954 

7782 



7121 

= 

0 










7122 

= 

0 










7 1 23 

= 

-0. 5651 

3209 

6413 

6430 

5307 

64 82 

3207 

0852 

• 

10 +1 

7124 

= 

0.4275 

0028 

7290 

4367 

7987 

3893 

2431 

0306 



7125 

= 

0 










7126 

= 

0 










7l 27 

= 

0.3029 

3416 

7269 

5682 

8108 

5679 

5437 

6506 



7l 28 

= 

-0.1028 

0329 

3795 

0334 

2611 

6141 

5109 

1571 

• 

10 +1 

7l 29 

= 

0. 5425 

4171 

2796 

6918 

2157 

8547 

6416 

2220 

• 

10 +1 

7l 210 

= 

0.1534 

0242 

6078 

6703 

1086 

6711 

9989 

5920 

• 

+ 

o 

71211- 

= 

-0.1576 

3473 

5858 

3826 

6589 

8937 

8096 

2028 

• 

10 

7130 

= 

0.3517 

3987 

5863 

0671 

3954 

7258 

1903 

7907 

• 

10” 1 

7131 

= 

0 










7132 

= 

0 










7133 

= 

0 










7134 


0 










7135 

= 

0 










7136 

= 

0 










7137 

= 

0. 6385 

8784 

3542 

5830 

8506 

8828 

9280 

0822 

• 

10 _1 

7138 

= 

0. 5086 

6724 

9055 

8144 

8754 

2910 

0714 

8603 

• 

10 -1 

7139 

= 

0.1770 

3179 

4727 

6675 

2427 

0314 

9422 

6269 



71310 

= 

0.1678 

1715 

6130 

4150 

9463 

9110 

6721 

5393 



71311 

= 

0.4538 

5629 

2579 

4244 

0722 

3753 

9295 

0401 

# 

10~ 2 

71312 

= 

0. 7129 

8936 

9976 

6658 

0243 

7127 

3010 

1115 

• 

10 -3 


i 


75 




TABLE 4. PATTERN FOR RKN-G-6(7)-10 
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TABLE 5. COEFFICIENTS FOR RKN-G-6( 7) -10 


O' 1 = 

0.1018 

5185 

1851 

8518 

5185 

1851 

8518 

5185 

G'2 = 

0.1527 

7777 

7777 

7777 

7777 

7777 

7777 

1778 

G 3 = 

0.2291 

6666 

6666 

6666 

6666 

6666 

6666 

6667 

g 4 = 

0.625 








O’ 5 = 

0.375 








0' 6 = 

0. 6666 

6666 

6666 

6666 

6666 

6666 

6666 

6667 

O' 1 = 

0.1666 

6666 

6666 

6666 

6666 

6666 

6666 

6667 

O' 8 = 

0. 9703 

7314 

8321 

7701 

1063 

1914 

4946 

5592 

0 9 = 

1 








O’ to = 

1 









010 

= 

0.1018 

5185 

1851 

8518 

5185 

1851 

8518 

5185 


020 

= 

0.3819 

4444 

4444 

4444 

4444 

4444 

4444 

4444 * 

10” 1 

021 

= 

0.1145 

8333 

3333 

3333 

3333 

3333 

3333 

3333 


030 

= 

0. 5729 

1666 

6666 

6666 

6666 

6666 

6666 

6667 • 

10 _1 

031 

= 

0 









032 

= 

0.1718 

75 








040 

= 

0. 8186 

9834 

7107 

4380 

1652 

8925 

6198 

3471 


041 

= 

0 









042 

= 

-0.3137 

9132 

2314 

0495 

8677 

6859 

5041 

3223 • 

10 +1 

0 43 

= 

0. 2944 

2148 

7603 

3057 

8512 

3966 

9421 

4876 • 

10 +1 

050 

= 

0 . 7840 

9090 

9090 

9090 

9090 

9090 

9090 

9091 • 

10 _1 

051 

= 

0 









052 

= 

0 









0 53 

= 

0.2906 

6985 

6459 

3301 

4354 

0669 

8564 

5933 


054 

- 

0.5921 

0526 

3157 

8947 

3684 

2105 

2631 

5789 • 

10" 2 

0 60 

= 

0. 8963 

7111 

8593 

3408 

1556 

3037 

7852 

6001 • 

10 _1 

061 

= 

0 









0 62 

- 

0 









063 

= 

0. 2041 

4673 

0462 

5199 

3620 

4146 

7304 

6252 


0 64 

= 

0.1424 

3014 

944 7 

6933 

0734 

2430 

1494 

4769 


065 

= 

0. 2304 

5267 

4897 

1193 

4156 

3786 

0082 

3045 


070 

= 

0.1018 

0041 

1522 

6337 

4485 

5967 

0781 

8930 


071 

= 

0 









072 

r= 

0 









073 

= 

0 









074 

— 

-0.3160 

4938 

2716 

0493 

8271 

6049 

3827 

1605 


075 

= 

0.1457 

9659 

0241 

0346 

8547 

9129 

9235 

7437 


076 

= 

0.2351 

1904 

7619 

0476 

1904 

7619 

0476 

1905 


080 

= 

-0.1674 

6873 

1305 

8845 

5834 

9712 

6386 

9845 • 

10 +1 

081 

= 

0 









082 

= 

0 












TABLE 5. (Continued) 


0. 7871 6193 2682 9395 6927 2722 4851 2526 • 10 +1 

0. 5392 4880 7631 6058 6316 9737 2744 6928 • 10 +1 

0.1544 8282 1 050 0 84 3 04 50 113 0 2 1 69 144 6 • 10 +1 

0.3255 5460 0653 6995 8589 4405 8687 3101 • 10 +1 

0.9924 5659 2893 1791 6591 1425 3844 6550 • 10 +1 

0.3470 3198 5621 8999 1428 6348 1503 3844 • 10 +1 

0 
0 

0.1593 8792 7828 8467 3788 8236 8550 9672 • 10 +2 

0.1140 4863 8019 8694 0109 1083 4461 7441 • 10 +2 

0.3469 8562 8685 7847 2962 3361 5811 0987 • 10 +1 

0. 7404 3819 3468 9437 2049 1200 7960 7603 • 10 +1 

0.1994 1980 7723 6203 6997 8978 8502 5319 • 10 +2 

0 . 634 9 3 71 3 6 9 80 1 9 6 7 41 73 4 3 8 8 5 7 84 4 7 7 9 • 10 _1 

0. 5371 0681 6961 8894 2565 7541 5995 7059 • 10" 1 

0 
0 
0 
0 

0.21 51 6847 1 682 4 6 71 944 5 4 911 4 7 3 9 2 647 

0.3420 7435 3923 4998 3829 7121 6887 6606 

0.2265 0729 0897 0719 7573 5894 5444 6475 

0.3047 2907 5218 4935 6793 3659 9050 3756 

0.1421 8987 3417 7215 1898 7341 7721 5190 


Yio 

— 

0. 5186 

8998 

6282 

5788 

7517 

1467 

7640 

6036 

• 

10 -2 

T 20 

5= 

0. 5835 

2623 

4567 

9012 

3456 

7901 

2345 

6790 

• 

10 -2 

y 21 

- 

0. 5835 

2623 

4567 

9012 

3456 

7901 

2345 

6790 

• 

10 -2 

730 

= 

0.1312 

9340 

2777 

7777 

7777 

7777 

7777 

7778 

• 

10' 1 

731 

= 

0 









10" 1 

732 

= 

0.1312 

9340 

2777 

7777 

7777 

7777 

7777 

7778 

• 

740 

= 

0.1775 

5681 

8181 

8181 

8181 

81 81 

8181 

8182 

# 

10 -1 

741 

= 

0 










742 

= 

0 










743 

= 

0.1775 

5681 

8181 

8181 

8181 

8181 

8181 

8182 


10 _1 

750 

= 

0.1925 

7489 

6694 

2148 

7603 

3057 

8512 

3967 

• 

751 

= 

0 









10 _1 

752 

= 

0.4028 

5268 

2912 

0077 

7831 

7938 

7457 

4623 

• 

753 

= 

0.1034 

9608 

5254 

4584 

6020 

0086 

9943 

4537 

• 

10 _1 

754 

= 

0.4201 

3351 

3931 

8885 

4489 

1640 

8668 

7307 

• 

10 -3 

7 60 

= 

0. 5015 

0891 

6323 

7311 

3854 

5953 

3607 

6818 

• 

10 _1 


083 = " 

084 = 

@85 = ~ 

08G = “ 
@87 = 
@90 ~ ~ 
@91 = 
@92 ~ 
@93 = “ 
@94 ~ 
@95 = - 
@96 = ~ 
@97 ~ 
@98 = ~ 
0100 ~ 
@101 = 
@102 = 
@103 ~ 
@104 = 
@105 ~ 
@106 = 
@107 = 
@108 ~ 
0109 = _ 


78 




TABLE 5. (Concluded) 


Yei = 0 

Y 6 2 = 0 


763 

= 

0.1283 

2080 

2005 

0125 

3132 

8320 

8020 

0501 



7 6 4 

= 

0.1427 

7669 

4823 

4784 

4920 

9443 

3614 

9015 

• 

10" 1 

7 6 5 


0.2947 

2859 

1024 

8873 

2118 

3617 

4799 

1378 

• 

10" 1 

770 

= 

0.1008 

4019 

2043 

8957 

4759 

9451 

3031 

5501 

• 

10 -1 

771 

= 

0 










772 

= 

0 










773 

= 

0 










774 

- 

-0.2032 

9218 

1069 

9588 

4773 

6625 

5144 

0329 

• 

10 -1 

775 

= 

0.8955 

5163 

6292 

3770 

3311 

7773 

8585 

1460 

t 

10" 2 

776 

= 

0.1517 

8571 

4285 

7142 

8571 

4285 

7142 

8571 

• 

10" 1 

Ygo 

= 

0. 5878 

8608 

9282 

4074 

8733 

2252 

3459 

2929 

• 

10 -1 

781 

= 

0 










782 

= 

-0. 6816 

7418 

1961 

6552 

8792 

8230 

9588 

5550 



783 

= 

0 










784 

= 

0 










785 

= 

0.3904 

9081 

2480 

6691 

0635 

4602 

0166 

2578 

• 

10 -1 

786 

— 

0.1320 

2858 

1017 

7748 

3299 

4568 

2086 

5647 



787 

= 

0. 9226 

1993 

4259 

5248 

2280 

2944 

4534 

3764 



790 

= 

0 . 84 3 5 

6399 

1262 

8663 

2857 

5687 

5609 

5560 

• 

10" 1 

791 

= 

0 










792 

= 

-0.1399 

94 3 7 

4195 

1129 

4662 

3845 

3271 

0268 

• 

10 +1 

793 

= 

0 










794 

= 

0 










795 

=: 

0 










796 

= 

0.1518 

2217 

0689 

4949 

2865 

2410 

0270 

4702 



797 

- 

0.1661 

2294 

3949 

3108 

2916 

6194 

1317 

3874 

• 

10 +1 

798 

= 

0.2535 

7326 

4223 

9613 

0665 

3317 

0496 

8283 

• 

10 -2 

Y 100 

= 

0. 5359 

5659 

84 20 

3665 

3630 

8565 

2940 

7254 

• 

10 _1 

7ioi 

= 

0 










7102 


0 










7l 03 

= 

0 










7104 

= 

0 










7i 05 

= 

0.1339 

3181 

3529 

5796 

4677 

4132 

9304 

7098 



7106 

= 

0.1144 

9729 

3051 

8163 

2233 

0537 

4835 

6721 



7l 07 

= 

0.1891 

5603 

3529 

7940 

7729 

1841 

4167 

3934 



7108 

=r 

0.7915 

0409 

1476 

6066 

6995 

7558 

1949 

0274 

• 

10 -2 

7109 

= 

0.9041 

5913 

2007 

2332 

7305 

6057 

8661 

844 5 

• 

10 -3 


79 





TABLE 6. PATTERN FOR RKN-G-5(6)-8 



8i 


I 














TABLE 7. COEFFICIENTS FOR RKN-G-5(6)-8 


cvi = 


15 ’ 


G' 2 ” 


5 ’ 


«3 = 5 ’ a 4 - 


10 


o _ ± , 

> O' 5 ~ A ’ q, 6 _ 7 » O- 7 - G 8 - 1 


010 

040 

053 

063 

070 

074 

080 

086 


_4_ 

15 


’ 020 ’ 021 " 1Q » 030 - 20 ’ 031 - 0 ’ 032 2 0 ’ 

9 . , 27 11 


40 


> 041 = 042 “ 0 ’ 043 " 77; ’ 050 - 4 g > 051 “ 052 “ 0 » 


5 _5 

> 054 4 g ’ 060 _ 


40 
27 112 


194 4 81 


» 061 062 0 » 


56 450 


64 827 
26 033 


* 064 


80 000 


194 481 


» 065 


24 544 


21 609 ’ 


41 796 
14 500 


236 575 

> 071 - 072 " 0 » 073 - - 38 313 


10 449 


> 075 _ 


275 936 


228 095 


45 279 


- 076 73 78g ’ 


81 ’ 081 “ 0 g2 0 83 ° ’ 0^ 3483 


250 . _ 160 

♦ 0 85 _ 7 


351 ’ 


2401 a - J_ 

> 0 87 in 


5590 ’ P87 10 * 


Yio 


81 729 _ 81 . 11 283 

740 ~ 640 ’ 741 “ ° ’ 742 3 2 00 ’ 743 1600 ’ 750 88 064 ’ 

3159 _ 7275_ _ _ _33_ _ 6250 

0 ’ 752 ~ 88 064 ’ 753 44 032 ’ 754 688 ’ 7g0 194 481 ’ 

3400 1696 6706 

7 62 -^63 — 0 » 7 64~ _ 194 4gl ’ Y65 " 64 83 7 ’ 7 70 4 5 279 

1 047 925 = 147 544 

771 - 772 - 773 - 0 ’ 774“ 4 940 997 ’ 7 75 " 19 g 209 ’ 


751 - 
761 = 


1 6 15 873 _ 31 _ _ _ „ , _ n 

776 1 874 886 ’ 7 80 3 0Q » 7 31 7 82 7 83 ’84 

_ _64_ .. 2401 _J_ 

785 585 ’ 786 7800 ’ 787 ~ 300 ‘ 




TABLE 8. APPLICATION OF THE VARIOUS FORMULAS TO PROBLEM I 



82 





















TABLE 9. APPLICATION OF THE VARIOUS FORMULAS TO PROBLEM II 



83 






























84 


NASA -Langley, 1974 


Figure 1. Orbit for the restricted problem of three bodies. 


